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Abstract 

A measurement of the inclusive ep scattering cross section is presented in the region of low 
momentum transfers, 0.2 GeV^ < < l2GeV^, and low Bjorken x, 5 ■ lO"'' < x < 0.02. 
The result is based on two data sets collected in dedicated runs by the HI Collaboration at 
HERA at beam energies of 27.6 GeV and 920 GeV for positrons and protons, respectively. 
A combination with data previously published by HI leads to a cross section measure- 
ment of a few percent accuracy. A kinematic reconstruction method exploiting radiative ep 
events extends the measurement to lower and larger x. The data are compared with the- 
oretical models which apply to the transition region from photoproduction to deep inelastic 
scattering. 
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1 Introduction 



Deep inelastic lepton-hadron scattering (DIS) is pivotal for the understanding of the structure 
of the nucleon and of the dynamics of parton interactions. Since the discovery of Bjorken scal- 
ing [1] and its violation [2] at fixed target experiments, DIS measurements have made essential 
contributions to in the development of the theory of strong interactions, Quantum Chromody- 
namics (QCD). Major progress in the exploration of strong interactions has been achieved at 
the electron '-proton collider HERA, operating at the energy frontier. Measurements performed 
at HERA are essential for predictions of the physics at the forthcoming proton-proton collider, 
the Large Hadron Collider (LHC). 

The high centre-of-mass energy of the ep scattering at HERA leads to a wide kinematic 
range extending to large values of the modulus of the four-momentum transfer squared, denoted 
2^, and to very small values of the Bjorken x variable. At the HERA beam energies of Eg = 
27.6 GeV for the electron and Ep = 920 GeV for the proton, Bjorken x values as small as 10"'' 
(10-^) are accessible for of lOGeV^ (0.1 GeV^). 

A salient feature of the structure function F2(x, Q^), discovered by the HI [3] and ZEUS [4] 
collaborations with the very first HERA data, is its strong rise for x ^ 0. In terms of parton 
distribution functions, this can be directly interpreted as a strong rise of the sea quark density 
towards small x. Similarly the increase of Fiix, Q^) with at fixed small x reveals a strongly 
rising behaviour of the gluon density towards low x. This is obtained in perturbative QCD 
(pQCD) analyses of DIS data [5-8] using the derivative dFildlwQ^, which is related to the 
gluon and quark densities as prescribed by the DGLAP evolution equations [9-13]. 

The DGLAP approach, in which only as In terms are summed, may not apply at lowest x 
values as terms involving powers of a, ln(l jx) become large. The parton dynamics at low x may 
be better approximated by different evolution equations, such as BFKL [14-16], CCFM [17-20] 
or non-linear equations [21-28]. The non-linear effects, arising due to the large gluon density 
and corresponding for example to gluon-gluon recombination, could tame the rise of F2 at low 
X. Further clarification of low x parton dynamics requires data of the highest precision, in a 
wide range of x and Q^. 

For < 2 GeV^, as the strong coupling constant a^iQ^) increases, the higher order cor- 
rections to the perturbative expansion become large and lead to the breakdown of the pQCD 
calculations. Measurements at low and low x thus probe this transition in which quarks and 
gluons cease to be relevant degrees of freedom. This onset of soft hadron physics is described 
by phenomenological, often QCD-inspired models (see [29] for a review). 

An attractive view of virtual photon-proton scattering has been developed with the colour 
dipole model [30]. It originated from the observation that in the proton rest frame, at low x 
the photon may fluctuate into a quark-antiquark pair with a lifetime oc 1 /x, long before the 
interaction with the proton [31,32]. Therefore the cross sections can be expressed as a product 
of the square of the wavefunction of the qq pair with a universal dipol-proton cross section. 
Another phenomenological model, used here, describes F2ix, Q^) based on the idea of self- 
similarity of the proton substructure at small x [33]. 

'Unless explicitly stated, the generic name "electron" is used throughout this paper to denote both electron and 
positron. 
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Access to the smallest x implies an extension of the measurements to high values of the 
inelasticity y where the cross section becomes sensitive to the longitudinal structure function 
Fl(x, Q^). This function completes the description of inclusive virtual photon-proton scattering, 
which involves transverse and longitudinal photon polarisation states. In the naive quark-parton 
model (QPM), Fi is zero, while in QCD it provides independent information [34] on the gluon 
distribution and may become correspondingly large at low x. 

This paper presents new measurements of the inclusive ep cross section in the range 0.2 < 
< 12GeV2 and 5 • 10"*' < x < 0.02. The data were collected with the HI detector in 
two e^p running periods with dedicated settings of the inclusive electron triggers. One data 
set (termed nominal vertex, "NVX") was collected in the year 1999 and corresponds to an 
integrated luminosity of 2.1 pb ^ The other was collected in the year 2000, with the interaction 
region shifted along the proton beam direction by 70 cm (termed shifted vertex, "SVX"), and 
corresponds to 505 nb '. 

Shifting the interaction region allows detection of the scattered electron at larger polar an- 
gles- which otherwise cannot be accessed in the main HI apparatus and thus provides accep- 
tance in the region <2 GeV^. In comparison to the previous HI measurement with a shifted 
vertex [35], an increased precision is reached using the higher luminosity of the new data and 
employing, in addition to the previous backward instrumentation of the HI detector, an up- 
graded Backward Silicon Tracker (BST). The vertex reconstruction using the electron track in 
the BST allows the kinematic range to be extended at low and low y. 

The measurement region is further extended towards lower and higher x values by ex- 
ploiting events with hard photons emitted coUinearly to the electron beam (Initial State Ra- 
diation or ISR). Such events are treated as ep events at an effectively reduced centre-of-mass 
energy. Unlike in the previous HI ISR analysis [36], the emitted photons are not explicitly 
detected, but the missing momentum is determined using momentum conservation. For this 
method the BST charged particle validation of the scattered electron is important to reduce 
the physics background from photoproduction events, in which the scattered electron escapes 
undetected in the electron beam direction. 

The measurement presented here is combined with previously published data [35,37] taken 
at Ep = 820 GeV in the region > 1.5 GeV^ (NVX97) and in the region > 0.35 GeV^ 
employing a shifted vertex technique (SVX95). The data sets are combined taking into account 
their systematic error correlations. The resulting accuracy reaches two percent precision in 
the bulk region of the measurement providing the most precise measurement in this kinematic 
domain. 

Data on F2 extending to low were published by the ZEUS Collaboration using a detector 
mounted near the beam pipe [38]. For > 2 GeV^, ZEUS data [39] from the 820 GeV 
operation of HERA are also available. 

The paper is organised as follows: In section 2 basic definitions are given. In section 3 
models are introduced which are subsequently compared to the data. In section 4 the methods 
to determine the DIS event kinematics and the principle of the cross section measurement are 

^In the HI coordinate system the z axis points along the outgoing proton beam direction termed forward di- 
rection. Therefore large electron polar angles 9,, close to 180° correspond to very small angles with respect to 
the incoming electron direction. The coordinate system is right-handed. The x (y) axis is directed horizontally 
(vertically). 
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presented. In sections the HI apparatus is described with emphasis on the components of key 
importance for the present measurement. Section 6 presents the event selection and reconstruc- 
tion, followed by section 7 on the Monte Carlo (MC) simulation of events. In section 8 a detailed 
account of the analysis techniques and uncertainties of the measurement is given, and the cross 
sections obtained from the 1999 and 2000 data are presented. In section 9 the data averaging 
method and the combination of the new data with the previous HI data taken at Ep = 820 GeV 
are presented. Section 10 is devoted to a phenomenological analysis of the x dependence of 
F2 and to extractions of the longitudinal structure function and in section 1 1 the data are 
compared to phenomenological models. A summary is given in section 12. 



2 Definitions 

In the low kinematic range of the present measurement, contributions from Z boson ex- 
change to neutral current deep inelastic scattering can be neglected. In the one-photon exchange 
approximation, the double differential cross section for neutral current DIS is given, in its re- 
duced form ar, by 

with the fine structure constant denoted a and f{y) = y^/[l -i- (1 - y)^]. The inelasticity y is 
related to Q^, x and the centre-of-mass energy squared, s = ^E^Ep, by 2/ = I sx. In the quark- 
parton model (QPM), x denotes the fraction of the proton momentum carried by the parton 
coupling to the exchanged boson. 

The DIS cross section, equation 1, is determined by two structure functions, F2 and F^. 
These are related to the cross sections for the scattering of longitudinally and transversely po- 
larised photons off protons, cti and ctt- At low x, the relationships 

(l-x)-cr^, (2) 



F2 = ^{\-x)-{ctl + (Tt), (3) 

hold to very good approximation. Positivity of the longitudinal and transverse scattering cross 
sections imposes the restriction < F^ < F2. Using the ratio R{x, Q^) 

R = ^ = -^, (4) 
ctt F2 - Fl 

the reduced cross section in equation 1 can be written as 



o-r = F2ix, <2') 



!-/(?/) 



R 



1 +R 



(5) 



For most of the kinematic domain, the reduced DIS neutral current scattering cross section is 
well approximated by the F2 structure function, since F^ leads to a sizeable effect only for large 
inelasticity values y. 
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The reduced cross section cr,. can be re-expressed as 



with the effective virtual photon-proton cross section 

o-f,p = aT^W-fm(TL. (7) 

The sum ctl -I- (jj is referred to as the total virtual photon-proton cross section, cr^°p, which 
is often expressed as a function of Q^- and of the invariant mass of the virtual photon-proton 
system, W . For small x, W can be calculated as = ^jQ^{l - x)/x, such that - sy. The 
total and the effective virtual photon-proton cross sections differ significantly only in the region 
ofhigh?/. 



3 Models 

The low X data presented here extend to low values of for which perturbative QCD is not 
applicable. The phenomenological models formulated for this transition region reproduce the 
W dependence of the y*p cross section, which is weak in the photoproduction region [40]. A 
steep increase towards large values of W develops in the perturbative region, which is equivalent 
to the rise of the proton structure function F2 towards low x at fixed Q^. 

In the context of the present measurement colour dipole models (e.g. [41-46]) are partic- 
ularly interesting because Fi and Fj = Fi - Fi are both described by a single characteristic 
dipole scattering cross section & combined with either the longitudinal or the transverse photon 
wavefunction. The squares of the wavefunctions of the qq fluctuations of longitudinally and 
transversally polarised photons are [30] 

6a ' 



' (8) 

3 

WAz, r, Q^) = -^Y^ e]{{\ - lz{\ - z))^K,{er) + m^^K^ier)] , 



i=\ 



where = m} + z{l - z)Q^, m,- (e,) is the mass (charge) of quark /, Kq{u) and ^i(m) = -BuKq 
are modified Bessel functions, r is the transverse separation of the qq pair and z denotes the 
fractional energy sharing between q and q. In this approach the cross sections ctt^l are obtained 
from integrals over the impact parameter space as 

(TL,Tix,Q^) = J d^rj^ dzWLAz,r,Q^)&(x,r^). (9) 



Colour dipole models differ by the chosen expressions for the cross section &. With the mea- 
surement extending into the region of high y one can confront the predictions of such models for 
the two structure functions with the data. As an illustration, the data are compared in this paper 
to two versions of the colour dipole model, the original version by Golec-Biemat and Wiisthoff 
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(GBW) [42] and a more recent model, based on the Colour Glass Condensate approach to the 
high parton density regime, by lancu, Itakura and Munier (IIM) [45]. 

Two further models are used in this paper in order to parameterise Fjix, Q^). The fractal 
model is based on the observation that the proton structure at low x exhibits self-similar prop- 
erties for different x and values. Two continuous, variable and correlated fractal dimensions 
are chosen to describe the self-similarity in x and [33]. In a more phenomenological ap- 
proach F2 is parameterised as x '^^^ \ These two models are also compared with the reduced 
cross section, cr^, after making assumptions on R. 



4 Measurement of the DIS Cross Section 



4.1 Reconstruction of Event Kinematics 



In the colliding beam experiments at HERA, the DIS event kinematics can be reconstructed 
using the measurements of the scattered lepton, the hadronic final state, or a combination of the 
two. This complementarity enlarges the kinematic coverage and provides an additional control 
of the systematic uncertainties. 

The energy of the scattered positron E'^ and its polar angle 6^ are used in the "electron 
method" to determine the kinematics via 



2£,-£;(l -cose,) 2£,-Z, 

= 2E, ^-IeT' ^'^^ 

, E'/ sin^ e, Ql 



Using energy-momentum conservation, the event kinematics can also be determined from 
the hadronic final state. An important quantity is the difference between the total energy and 
the total longitudinal momentum 

E-P, = K{\- cos e,) + J] (Ei - P,j) = + 2, , (12) 

i 

where E, (P^j) is the reconstructed energy (longitudinal component of the momentum) of a 
particle i from the hadronic final state. In the reconstruction masses are neglected for both the 
positron and the hadronic final state particles. The measured E-P^ is insensitive to losses in the 
proton beam direction and is thus only weakly affected by the incomplete reconstruction of the 
proton remnant. For non-radiative events, the relation E-P^ ^ lEg holds. This allows lE^ - Sg 
in equation 10 to be replaced by l^i, and leads to the introduction of the variable [47] 

yh = ^. (13) 



For events in which a photon is emitted coUinearly to the incoming positron, the radiated 
photon is not reconstructed in the sub-detectors used to calculate E-P^. In this case {E-P-)/2 
is equal to an "effective" incident positron beam energy, reduced relatively to the nominal beam 
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energy by the momentum carried by the radiated photon. This is employed in the 2 method, for 
which lEe in equation 13 is substituted by the measured E-P^ [48] 



(14) 

For this method, is calculated by replacing ye in equation 1 1 by j/s, and Bjorken x is calcu- 
lated by substituting i/g, Ql and lEg by i/i, Q\ and E-P,, respectively'' 

By using a consistent set of the variables x^,, yi. and Q^, the measurement also correctly recon- 
structs the kinematics for events with initial state QED radiation. Therefore, the method covers 
lower and higher x values, which become accessible due to the reduced centre-of-mass en- 
ergy for these events. 

The total transverse momentum of the hadronic final state is 



pii _ 



(16) 



where Ft,/ is the transverse momentum vector of the particle i and the sum runs over all parti- 
cles. Pj is rather insensitive to particle losses coUinear to the beam for a wide range of y. The 
combination of Pj and S/, defines the hadronic scattering angle 

Oh 2/i 

tan — = — , (17) 
2 

which, within the QPM, follows the direction of the struck quark. 

In this analysis, both the electron and the 2 methods are used for the cross section measure- 
ment. The electron method provides the better resolution in x for large inelasticities y > 0.1, 
but the resolution degrades as l/y. Use of the S method extends the measurement down to 
y ~ 0.002. Below this y value, losses along the proton beam direction become important and 
are difficult to estimate. The S method as is used here noticeably increases the kinematic cov- 
erage towards low and high Bjorken x due to initial state QED radiation. 



4.2 Determination of the DIS Cross Section 

The measurement of the double differential cross section is performed in bins of x and Q^, or 
y and Q^, depending on the region in the kinematic phase space, as shown in figure 1. The bin 
sizes and shapes as well as methods used for the kinematic reconstruction are chosen based on 
the following prescription: 

• In Q^, a binning equidistant in logjo is chosen with eight bins per decade, as in previous 
HI publications [37]. This binning reflects the good resolution of the HI detector. 

^Note that in previous HI publications the nominal positron beam energy was used instead of {E-P-)/2 in the 
calculation of .yj;. The method of x calculation used here is called the IE method in [48]. 
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• The X and values at which the measured double differential cross section is quoted, 
also referred to as bin centres, are placed at an approximately logarithmic average value 
within the bin boundaries for the x and binning, and at the linear average for the y 
binning. 

• For high y > 0.6, the electron method has an excellent kinematic resolution. In this 
region, the measured cross section is sensitive to the longitudinal structure function Fl, 
which leads to a rapidly changing shape of the reduced cross section as a function of y. 
Therefore a fine binning, linear in y, is chosen for y > 0.6: two y bins are used for each 

interval with boundaries at y = 0.85, 0.75 and 0.6. 

• For y < 0.6 the binning is defined in x. The default x binning is equidistant in logjpX 
with five bins per decade, as chosen previously [37]. The transition between the x and y 
binning is defined by the y value of the nominal bin centre, y^, for the transition bins: for 
yc > 0.6, the bin is combined with the nearest y bin and for y^ < 0.6 it is combined with 
the nearest x bin. 

The resolution in each bin is checked using a Monte Carlo simulation. Two variables are cal- 
culated for this purpose, the purity P = A^recgen/A^rec and the stability S = A^recgen/A^gen, where 
A'rec (A'^gen) IS the total numbcr of reconstructed (generated) Monte Carlo events in the bin and 
A'^recgen IS the number of events which are both generated and reconstructed in the same bin. The 
purity and stability are calculated for both the electron and the E methods. For the cross section 
measurement the method with the higher purity is used. The choices are illustrated in figure 1 . 
The purity and stability typically exceed 50%. If either the purity or the stability is below 25% 
in a bin for the chosen reconstruction method, the bin is combined with an adjacent bin. Bins 
with larger sizes can thus be created at the acceptance edges as shown in figure 1. 

The 6e = 176.5° and Og = 178° lines in figure 1 indicate the approximate angular acceptance 
limits of the HI detector for the nominal and the shifted vertex positions, respectively. In each 
plot measurement bins below lines are visible. The measurement in these bins becomes 
possible using the 2 method for events with initial state photon radiation which effectively 
reduces the centre-of-mass energy. These bins are further referred to as ISR bins. The E method 
cannot be used at high y, where its resolution is poor, leading to large migrations of nominal 
energy events into the ISR bins and thus to purities below the accepted value. This causes the 
gap between the ISR and electron method bins at high y. 

The calculation of the reduced double differential ep cross section is performed by cor- 
recting the data using the MC simulations. The following formula is applied to each analysis 
bin 



Here, {Xc, Q^) is the bin centre, A^^ata is the number of data events, A^bg is the number of back- 
ground events, estimated using MC simulations, A and e are the detector acceptance and ef- 
ficiency, Xdata is the integrated luminosity, 6^ are QED radiative corrections, and Cbc are the 
corrections for finite bin size effects. The radiative and bin centre corrections can be deter- 
mined using the Monte Carlo simulation. In this case, equation 18 becomes 




(18) 




(19) 
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where A^^mc is the number of signal MC events and Xmc = A^gen/o"gen is the Monte Carlo lu- 
minosity. Here A^gen denotes the total number of generated MC events and cTgen is the total 
integrated cross section for the MC generation. The quantity cr^'~^ (x^, Qfj is the reduced double 
differential cross section at the bin centre calculated at the Born level with the same structure 
functions as are used in the MC generation. 

The correction for the detector acceptance using Monte Carlo modelling requires the cross 
section model used in the simulation to be sufficiently close to the data, such that migrations be- 
tween the bins are well reproduced. The cross section model should also describe the kinematic 
region outside the measurement range, in particular at low y and low Q^, to account for radiative 
corrections and long range migrations. In practice, this is achieved using an iterative MC event 
reweighting procedure which converges after one iteration for the measurement region. First, 
the double differential cross section is measured following equation 19 using an initial approx- 
imation for the MC input cross section. Next, the measured double differential cross section is 
fitted with a new parameterisation using the fractal model and the analysis of the Monte Carlo 
events is repeated with an additional weight factor, equal to the ratio of the new to the initial 
double differential cross sections in each simulated event. For the reweighting, the event kine- 
matics are calculated using the generated x and variables at the hadronic vertex, such that 
corrections due to radiation from the lepton line are properly accounted for. This reweighting 
procedure is used for the measurement region. For the high x > 0.02 domain, which lies outside 
the measurement region, the ALLM parameterisation [49] is used. 



5 HI Detector 

5.1 Overview 

A complete description of the HI detector is given in [50,51]. Here the components used for 
the present measurement are discussed. In section 5.2 the detectors for the scattered electron 
measurement are described in detail. A schematic view of the HI detector is given in figure 2, 
in which a typical low event is shown. 

Around the interaction region a set of tracking chambers, surrounded by electromagnetic 
and hadronic calorimeters, operates in a solenoidal magnetic field of 1.16 T. The tracking sys- 
tem is subdivided into forward, central and backward tracking devices. The nominal interaction 
point of the electron and proton beams lies about in the middle of the Central Tracker, at the 
origin of the coordinate system. The interaction vertex positions have an approximately Gaus- 
sian distribution in z with cr. « 10 cm. The calorimetry system consists of the Liquid Argon 
calorimeter (LAr) covering the central and forward directions and the lead-scintillator spaghetti 
calorimeter (SpaCal) [52-55] measuring particles scattered backwards. 

The Central Tracker consists of four drift chambers, two multi-wire proportional chambers 
(MWPCs) and a silicon tracking device. The largest tracker components are the two concentric 
drift chambers, CJCl and CJC2, which have sense wires strung parallel to the beam axis with the 
drift cells inclined at about 30° with respect to the radial direction, such that the drift direction 
of ionisation electrons is approximately perpendicular to the wire plane. The charge deposits 
are read out from both ends of each wire, providing particle identification via ionisation energy 
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loss and an approximate determination of the z coordinate via the charge asymmetry between 
the two wire end signals ("charge division"). 

Tracks found in the CJC are linked to the hits found in two chambers equipped with wires 
strung around the beam axis, following polygonal support structures, dedicated to the precise 
measurement of z coordinates. The inner z chamber (CIZ) is located inside CJCl and the outer 
z chamber (COZ) lies between CJCl and CJC2. To reduce the number of acceptable combi- 
nations with the CJC, the z chambers also determine a (p coordinate using the charge division 
measurement. The tracks are further constrained by linking to hits in the central silicon tracker 
(CST) [56]. The CST consists of two layers of double-sided silicon strip detectors surrounding 
the beam pipe, covering an angular range of 30° < 6 < 150° for tracks passing through both 
layers. 

The two cylindrical proportional chambers, the CIP mounted inside CIZ, and the COP lo- 
cated between the COZ and CJC2, are used together to identify tracks pointing to the interaction 
vertex and thus to reduce background at the trigger level. A combined CIP-COP signal is used 
in coincidence with the SpaCal to trigger events with low E'^ (see section 6.1). 

The LAr calorimeter [57], mounted in a large cryostat, is used in this analysis for the mea- 
surement of the hadronic energy. The angular coverage of the calorimeter is 4° < < 154° for 
an interaction vertex at z = 0. The calorimeter consists of an electromagnetic section with lead 
absorbers (20-30 radiation lengths) and a hadronic section with steel absorbers. The total depth 
is between 4.5 and 8 hadronic interaction lengths. The LAr calorimeter is divided along the z 
direction into wheels. The electromagnetic section has eight wheels while the hadronic sec- 
tion has seven. The calorimeter has a high degree of spatial segmentation with a total of about 
45000 cells. Its hadronic energy resolution, as determined in test beam measurements [58], is 
aslE ^ 50%/ V£/GeVe2%. 

Two electromagnetic crystal calorimeters, a photon tagger (PT) and an electron tagger (ET), 
located at z = -103.1 m and z = -33 m, respectively, are used to monitor the luminosity via 
the measurement of the Bethe-Heitler process ep yep. The luminosity corresponding to 
the main interaction region can be separated from the additional ("satellite") interaction regions 
using information from the scintillator hodoscopes of the time-of-flight system (TOP) and from 
the HERA proton pick-up (PPU) monitor, a 34 cm long stripline device located at -3 m from 
the interaction point. The ET can be used to measure the scattered electron in photoproduction 
processes, with < 10"^ GeV^ and Q.2 <y < 0.7. The PT detects photons radiated coUinearly 
to the incoming electron direction. 

5.2 Backward Detectors 

The measurement of the inelastic ep scattering cross section at low relies on the identification 
of the scattered electron in the backward part of the HI apparatus. The energy of the scattered 
electron is measured in the SpaCal calorimeter. For the low region under study, 6^ lies 
outside the angular acceptance of the Central Tracker. The polar angle of the scattered electron 
can, however, be measured either by the Backward Silicon Tracker (EST), based solely on the 
electron track, or by a combination of the less precise Backward Drift Chamber (BDC) signal 
with the hadronic final state vertex, as reconstructed using the Central Tracker. The redundancy 
of the angular measurements provides additional cross checks over a large angular range, whilst 
the BDC extends the polar angle coverage to larger 6e. 
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5.2.1 Backward Silicon Tracker 

The BST in the configuration installed in 1999 [59] is schematically shown in figure 3. It con- 
sists of eight planes (disks) and 16 azimuthal sectors. The planes are mounted perpendicularly 
to the beam axis and are arranged in two modules, BSTl and BST2, of four planes each. A first 
version of the BST with four planes is described in [60]. 

Each BST plane is equipped with 16 wedge shaped, single sided, double metal, silicon strip 
sensors of 250 fxm thickness. Each sensor contains 640 sensitive p strips which are concentric 
around the beam axis with a pitch of 96 /im. The signals are amplified and temporarily stored by 
five on-detector front-end chips, called Analogue Pipeline Chips [61] (APCs), until a readout 
instruction is received [62]. Using these "r sensors" (figure 4a) the track polar angle can be 
determined. The acceptance range of the BST for the nominal vertex position is 164° < 6e < 
176°. 

In addition to the r sensors, each plane contains one single sided, single metal, silicon 
strip sensor, in the azimuthal sector 45° < (p < 67.5° mounted behind the r sensor. This "m 
sensor" has 640 sensitive strips parallel to the reference edge of the sensor with a pitch of 75 jum 
(figure 4b). It thus measures hits in u coordinate space defined by m = r sin0„, where 0„ is the 
azimuthal angle with respect to the reference edge of the sensor. Combining the information 
from r and u sensors, it is possible to measure the transverse momentum and determine the 
charge of a track in the BST. This feature is used in this analysis to cross check the simulation 
of photoproduction background. 

During data taking an online hit finding is performed. This takes into account individual 
pedestals of each channel, which are dynamically updated. Coherent shifts in the amplitude of 
groups of strips, so called "common mode", are also corrected for. For reconstructed tracks, the 
most probable signal-to-noise values for the hits is about 15 for the r sensors and 30 for the u 
sensors. The single hit resolution is 20 (15) ijm for the r (u) coordinate. 

5.2.2 SpaCal and BDC 

The SpaCal calorimeter covers the polar angle range of 153° < 6 < 177° as measured from the 
nominal z vertex position. It consists of an electromagnetic section [53-55] with 1192 cells of 
size 4.05 x 4.05 x 25 cm^ in front of a hadronic section with 136 cells of size 11.9x11.9x25 cm^^ . 
The total amount of passive material traversed by particles from the interaction vertex up to the 
SpaCal is of the order of one radiation length. The electromagnetic section comprises 27.5 radi- 
ation lengths and provides an electromagnetic energy resolution of cte/E = 7%/ V£7GeV©l%. 
The hadronic section [52] is used for a coarse hadronic energy measurement and to distinguish 
hadronic from electromagnetic showers. The whole calorimeter comprises 2 hadronic interac- 
tion lengths. The energy resolution for hadrons amounts to cte/E ~ 60%/ ^/E/GeV. 

The SpaCal cells consist of lead sheets with embedded scintillating fibres. The fibres from 
each cell are bundled together and attached via light mixers to photomultiplier tubes (PMTs). 
The stability of the PMT gain can be checked using a dedicated LED system. 

The backward drift chamber [63] is mounted in front of the SpaCal and has the same angular 
acceptance. It consists of four double layers, each of them divided azimuthally into eight sec- 
tors. A three dimensional view of a section of the BDC is given in figure 5. The sense wires are 
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Table 1: Overview of the main trigger conditions used for the NVX and SVX analyses, r^p is 
the radial coordinate of the SpaCal cluster. 

strung perpendicularly to the beam axis and are fixed at the sector edges leading to an octagonal 
geometry with almost radial drift directions. The drift cells are 1 cm wide in the inner region 
and 3 cm wide in the outer. At the transition from the inner to the outer region a special cell 
is introduced with 0.5 cm drift distance at the inner side and 1.5 cm drift distance at the outer. 
The cells within one double layer are shifted by half a drift cell in the radial direction to solve 
the inner-outer hit assignment ambiguity. The double layers are rotated by 1 1.25° with respect 
to each other to allow for reconstruction of the azimuthal coordinate. In addition, this reduces 
the efficiency losses at the sector edges. The radial resolution for minimum ionising particles is 
400 //m. The resolution in the azimuthal direction is about 2 mm. 



6 Data Collection and Reconstruction 

The HI detector uses a multi-level trigger system for data collection in which two hardware 
trigger levels are followed by a software filter farm. After improvements of the detector cali- 
bration and the reconstruction code, the data are reprocessed offline. This section describes the 
ffrst analysis stages, including the online data selection and the reconstruction algorithms. 

6.1 Online Event Selection 

The online trigger conditions used in this analysis (table 1) are based primarily on a localised en- 
ergy deposition in the electromagnetic section of the SpaCal (inclusive electron trigger). Three 
difi'erent energy thresholds are used. The trigger condition based on the lowest energy threshold 
(S9) collects events at the highest y. Since a significant background contamination is present at 
low energies and radii, the inner SpaCal region is excluded from S9. To maintain an acceptable 
trigger rate, S9 contains an extra condition requiring the pattern of hit pads in the CIP-COP 
proportional chambers to be consistent with at least one vertex pointing track. This condition is 
similar to the requirement of a reconstructed track from the hadronic final state in the Central 
Tracker. 

6.2 Track and Vertex Reconstruction in the Central Tracker 

The track reconstruction in the Central Tracker is initiated from the CJC hit measurements. 
Initially, candidate trajectories are found in the xy plane using a fast circle fit algorithm [64]. 
The z coordinate is added to the tracks based on charge division information. A linear fit in 
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S - z space is perfomied where the S coordinate measures an arc length of the track in the xy 
projection. Next, the tracks are fitted to a common vertex in the xy plane. At this stage detailed 
corrections are applied for multiple scattering in the detector material and for magnetic field 
non-uniformity. For the tracks which are determined by the fit to originate from a common xy 
vertex, a combined S - z fit is performed to determine an initial approximation of the vertex z 
position and of the track polar angles. 

The vertex-fitted CJC tracks are then combined with hits found in the z chambers, employing 
a robust global minimisation technique [65]. This link improves the z vertex resolution from 
about 1 cm to 1 mm. The track 9 resolution is correspondingly improved from approximately 
lOmrad to 1 mrad. For the NVX sample, where the z coordinate of the interaction vertex is 
around zero, the CJC tracks are also combined with hits found in the CST resulting in a vertex 
resolution of about 0.1 mm. CST hits are not used for the SVX sample since the interaction 
vertex position is outside the CST acceptance. 

6.3 Reconstruction of the Electron Kinematics 

6.3.1 Energy Reconstruction in the SpaCal Calorimeter 

The reconstruction of the scattered electron kinematics is based on the measurement of a deposi- 
tion of energy, termed a cluster, found in the electromagnetic SpaCal. The clustering algorithm 
searches for a cell with a local maximum in energy. The cluster is then built around this seed cell 
by adding neighbouring cells with energies above the noise threshold. The centre-of-gravity of 
the cluster is determined based on all associated cells using a logarithmic energy weighting. To 
suppress background from hadrons and from decays of yy with the photons reconstructed 
in a single cluster, a cluster radius estimator, 7?iog, is used based on logarithmic energy weight- 
ing. The background from hadronic particles is further suppressed using a cut on the energy 
deposit, Ehad^ in the hadronic section of the SpaCal behind the electromagnetic cluster. 

The electron candidate cluster is required to be associated to a track in one of the backward 
trackers, BDC or BST, in order to reduce background from backward photons and to measure 
the polar angle Qe accurately. For the determination of the trajectory of the scattered electron 
is assumed to be a straight line in r - z coordinate space. 

6.3.2 Track Reconstruction in the BDC 

The BDC reconstruction of the electron scattering angle Qe starts from the line connecting the 
SpaCal cluster and the Central Tracker vertex as an initial approximation. The scattered electron 
azimuthal angle 0e is taken from the SpaCal cluster centre-of-gravity. Only the BDC hits in the 
octant containing 0e are used for the reconstruction. 

The determination follows from a minimisation procedure. A least squares track fit com- 
bines the Central Tracker vertex, the SpaCal cluster centre-of-gravity energy, and all BDC mea- 
surements in a corridor of variable size Ar around the current best estimate of the track direction. 
Initially, the corridor has a size of 5 cm. It is gradually reduced with improved track parameters 
to about five times the BDC resolution. The SpaCal cluster is considered to be linked to the 
BDC track segment if there are at least four hits from the eight layers remaining at the final 
iteration and if the radial distance between the track projected to the SpaCal z coordinate and 
the SpaCal cluster is less than 2.5 cm. 
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6.3.3 Track Reconstruction in the BST 



The reconstruction of the electron track in the BST uses the azimuthal location of the SpaCal 
cluster. The three adjacent sectors which in azimuth are closest to 0e are selected. The r co- 
ordinates of all BST hits in the selected sectors are projected along the line defined by the hit 
and the SpaCal cluster to the BST plane closest to the SpaCal. A clustering of the projected 
hits in this plane is then performed using a histogram technique. The line connecting the po- 
sition corresponding to the peak in the histogram and the SpaCal cluster is used as an initial 
approximation for the track. 

The track finding then proceeds using an iterative minimisation technique with robust rejec- 
tion of outliers, similar to the BDC reconstruction. All hits in the selected sectors are included 
into a least squares minimisation. The contribution of each hit is weighted with an exponential 
suppression factor, which depends on the distance from the hit to the track, and on an additional 
parameter, which defines the width of an effective corridor around the track. For the first iter- 
ation, the width of the corridor is equal to the SpaCal spatial resolution. For further iterations 
the width is gradually reduced until it reaches five times the BST spatial resolution. The event 
vertex z coordinate is given by the distance of closest approach of the BST track to the beam 
line. 

For the sector equipped with the u strip detectors, the reconstruction of the azimuthal coor- 
dinate is also performed. At least three u hits associated to linked r hits are required. If multiple 
u hits per plane are found, all possible track combinations are formed and the one best match- 
ing the SpaCal cluster is selected. To determine the space points, the u hits are combined with 
the r hits extrapolated along the r track to the z position of the u sensor. Then the transforma- 
tion (r, u) => (x, y) is performed. A circle fit including the position of the interaction vertex in 
{x, y) determined by the beam spot size of 150 yum in x and 60 jjm in y, yields the curvature and 
therefore charge and the transverse momentum of the particle. 

6.4 Reconstruction of the Hadronic Final State 

The reconstruction of the hadronic final state uses information from the central tracker and the 
LAr and SpaCal calorimeters, excluding a cone in the SpaCal calorimeter around the electron 
candidate cluster. The cone axis is defined by the vertex position and the centre-of-gravity of the 
SpaCal cluster. The cone radius is 20 cm at the surface of the SpaCal electromagnetic section. 
The energy of the cells inside the cone is excluded from the hadronic final state calculation for 
both the electromagnetic and hadronic sections of the SpaCal. 

Tracks pointing to the backward part of the HI detector are excluded from the hadronic final 
state. Instead, the reconstructed SpaCal clusters outside the electron isolation cone are used. In 
the central region, the Central Tracker and LAr signals are linked for each particle by matching 
the measurements in each detector. For energies below 2 GeV, the tracker information is used 
while for higher energies the calorimeter information is used, as it provides the better energy 
resolution. 

The determination of S/, is affected by the presence of extra activity in the calorimeters. The 
bias is particularly strong for small S/, and thus small i//,. For the SpaCal, this extra activity can 
be induced by the scattered electron, with some energy leaking outside the isolation cone or by 
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a radiated photon emitted at a large angle. The contribution of these sources of extra activity to 
E/j is proportional to 2^ to good approximation. To reduce the influence of these effects, 10% of 
is subtracted from the total SpaCal E/j. If the result is negative, the SpaCal contribution is set 
to zero. This procedure reduces the contribution of SpaCal to E/, to a negligible level for low y 
events, as is expected from the event kinematics. 

Channels afl'ected by electronic noise in the LAr are identified event by event using a ded- 
icated topological algorithm. LAr cells with an energy below 0.4 GeV (0.8 GeV), which are 
separated from other cells by more than 40 cm (20 cm) in the central (forward) region of the 
calorimeter are classified as noise and excluded from the and P\ calculations. 

7 Monte Carlo Simulations 

In the simulation, DIS events are generated using the DJANGOH 1 .4 [66] event generator which 
includes leading order QED radiative efi'ects as implemented in HERACLES [67]. For the event 
generation, leading order parton distribution functions define F2 while Fi is set to zero. The 
structure functions are subsequently reweighted to the fractal model parameterisation of F2 and 
to Fl following the procedure described in section 4.2. The final state parton showers are sim- 
ulated using the Colour Dipole Model [68,69] as implemented in ARIADNE 4.1 [70]. Events 
with a very low mass of the hadronic final state {W < 5 GeV) are simulated using SOPHIA [71], 
which includes a detailed description of low mass final states, including the resonance region. 
The fragmentation into hadrons is performed with JETSET 7.4 [72]. Photoproduction back- 
ground is generated with the PHOJET 1.6 [73,74] program, which uses a two-component dual 
parton model [75] including difFractive processes and vector meson production. 

The simulation of QED radiative corrections includes photon emission from the lepton. Ra- 
diation from quarks, which is estimated to be small for low x, is not simulated. The simulation 
of QED radiative corrections is checked using the analytical calculation package HECTOR [76]. 
An agreement to better than 0.5% is found in the kinematic range of this measurement. 

The generated events are passed through a simulation of the HI detector response based 
on the GEANT3 [77] package. Tracing of the particles in the trackers up to the calorimeters 
is based on a detailed description of the detector material. The response of the calorimeters 
to electromagnetic particles is simulated using a fast shower parameterisation technique [78], 
while the hadronic response is simulated using GHEISHA [79]. 

The level of noise and beam related background in the calorimeters is determined using 
events from dedicated runs with random triggers which are overlaid on the simulated events. 
Spurious hits in the BST are added to the simulation based on randomly triggered events. 

The MC events are subjected to the same reconstruction and analysis procedure as the data. 
Also, for consistency of the analysis, the calibrations of the SpaCal and the LAr, as well as the 
BST and BDC alignments, are performed for the reconstructed MC events in the same way as 
for the data. 
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Analysis 


Trigger 


Description 


NVX(zvt, ~Ocm) 


NVX-BST 


SO 


Main analysis 


NVX-SQ 








SVX (zvtx 


~ 70 cm) 


SVX-BST 


SO 


Main analysis 


SVX-BDC 


so, S3 


Extension to larger 



Table 2: Inclusive analyses of DIS data. The trigger conditions used to collect the data sets are 
described in section 6.1, table 1. 

8 Data Analysis 

At low the DIS cross section is large, and for the available integrated luminosity for this 
analysis the statistical uncertainty of the measurement becomes smaller than the systematic un- 
certainty. For low inelasticities, corresponding to a large fraction of the measured phase space, 
the scattered electron energy is large and background contributions are negligible compared to 
the genuine DIS signal. In this region a set of selection criteria is imposed which is sufficient 
to reconstruct the event kinematics in the least biased way. Whenever possible the electron 
trajectory is reconstructed using the EST alone and only the SpaCal is used for triggering. 

Events with the scattered electron outside the EST acceptance are reconstructed using the 
EDC and the Central Tracker vertex. The analysis is also extended to the highest accessible y 
values for which the precision is limited by the uncertainty of photoproduction background. In 
this region several additional electron identification criteria are imposed in order to minimise 
the systematic uncertainty. 

For the two data samples, NVX and SVX, a total of four separate analyses is performed 
as summarised in table 2. The analyses diff"er in the triggers and in the method employed for 
reconstructing Qg. The main kinematic region of the NVX-EST data set, with Qe measured in 
the EST, is analysed based on the trigger SO. An extension to 0.75 <i) < 0.85 is achieved using 
the trigger S9 (NVX-S9) and requiring signals in both tracking detectors, EDC and EST. For 
the SVX sample, the main region of the phase space is covered by a EST-based analysis with 
the trigger SO (SVX-EST). An extension to Be = 178° is achieved by adding data collected with 
the trigger S3 and including events with 6^ measured by a combination of the Central Tracker 
vertex and EDC information (SVX-EDC). 

The measurement is verified by performing a number of cross check analyses exploiting the 
redundancy in the kinematic reconstruction and the large overlap of the kinematic regions of 
diff"erent data sets. The data reconstructed with the EST are compared with those reconstructed 
with the EDC. The results of the electron method are cross checked with those of the 2 method. 
Moreover, the measurement based on the shifted vertex sample is compared to that based on the 
nominal vertex sample. 

In the following a detailed description of the different analyses is given. Further information 
can be found in [80-83]. 
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Description 



Cut 



Common cuts 



Scattered electron energy 


E[>1 GeV; 4 GeV (NVX-S9) 


Vertex z position 


kvtx - ZnomI < 35 Cm 


SpaCal cluster radius 


i?iog < 4 cm 


Hadronic energy fraction 


£h,d/£:. <o.i5 


Pt balance 


P^^/P^ > 0.3 


Electron method cuts 


E-P, balance 


E-P, > 35 GeV 


BST analysis cuts for NVX-S9, NVX-BST and SVX-BST 


EST validation 


Mink BST > 2 (NVX); 3 (SVX) 


BST-SpaCal radial match 


lAresT-spaCail < 1.5 cm 


BST noise 


Mhit total < 120 (NVX); 200 (SVX) 


BDC analysis cuts for SVX-BDC 


BDC validation 


Mink BDC ^ 4 


BDC-SpaCal radial match 


lArBDc-SpaCail < 2.5 cm 


Central Tracker vertex 


Mtrack ^ 1 




> 0.03 


Additional NVX-S9 analysis cuts 


BST-BDC radial match 


|ArBST-BDcl < 0.75 cm 


BST-CT zvtx match 


kvTX.BST - ZVTX,CT /C < 5.0 


Central Tracker vertex 


Mtrack ^ 2 



Table 3: Selection criteria used in the analysis. 



8.1 Event Selection 
8.1.1 Criteria 

An overview of the selection criteria used in the different analyses is given in table 3. The 
background from non-ep interactions is suppressed by requiring the event vertex (zvtx) to be 
reconstructed within a distance of ±35 cm from the average z position (znom)- In order to be 
identified with the scattered electron, the highest energy cluster in the electromagnetic SpaCal 
section"^ has to satisfy the following criteria: (i) the cluster centre-of-gravity lies in the region 
of high efficiency of the corresponding trigger; (ii) the transverse cluster radius is consistent 
with an electromagnetic particle, R\og < 4 cm; (iii) the energy deposition in the hadronic SpaCal 
section behind the cluster is small, Ehad/^^^ < 0.15; (iv) depending on the analysis, the cluster is 
validated by a BST or a BDC track segment. If the highest energy cluster does not satisfy one 
of these cuts, the next highest energy cluster is used. This procedure is repeated for up to three 
clusters with energies above 7 GeV, or 4 GeV (NVX-S9). 

The further event selection is based on a global balance between the hadronic final state 
and the electron. Events for which the hadronic final state is poorly reconstructed are rejected 

"^For the S9 analysis (table 2) the cluster with the maximum transverse momentum is chosen instead of the 
highest energy cluster 
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by demanding that the total hadronic transverse momentum Pt^ be at least 30% of the electron 
transverse momentum Pj. This efficiently removes migrations from very low y, which lie out- 
side the measurement region. Events with large initial state radiation are excluded from the 
electron method measurement by requiring E-P- > 35 GeV. This condition is not used for the 
E method, which takes QED radiation properly into account. 

The BST analyses include requirements on the minimum number of BST hits linked to the 
electron track (Mink bst) and on the matching of the BST track extrapolated to the z position of 
the SpaCal cluster, ArBsx-SpaCai- Similarly, for the BDC based analyses, a minimum number of 
linked BDC hits (Mink bdc) and radial BDC-SpaCal matching (ArBoc-spaCai) are required. In 
addition, the BST analyses require a low level of noise by cutting on the variable Nhn totab the 
total number of BST hits. The BDC analyses demand the presence of at least one central track 

(Mtrack)- 

The S9 analysis extends the measurement to low E'^, corresponding to high y, where the 
largest uncertainty stems from the large photoproduction background. To suppress this back- 
ground, both the BST and BDC track segments are required to pass the same criteria as in 
the other analyses. In addition, a tight matching condition is applied for the two trackers using 
A^BST BDC, the radial distance between the BDC and BST tracks calculated at the BDC plane, as 
well as kvTx,BST - Zvtx.ctI/o', the distance in z between the BST vertex and the Central Tracker 
vertex position divided by the uncertainty of this difference. Finally, to ensure a high trigger 
efficiency for the analysed sample, at least two central tracks must be reconstructed. 

8.1.2 Efficiency Determination 

The efficiencies of the triggers are determined using independently triggered data samples. For 
the SpaCal trigger conditions, most of the cells show a high (> 99.5%) efficiency above the 
rather sharply defined threshold, see figure 6. A few cells are identified which show high thresh- 
olds. They are excluded from the analysis by applying geometrical cuts on the electron impact 
point reconstructed at the calorimeter surface, which is calculated using 6^ and 0g. The effi- 
ciencies of the CIP-COP conditions employed in the S9 trigger (see section 6.1) are studied 
as functions of E'^ and the track multiplicity. Since the average reconstructed track multiplicity 
increases with Q^, the inefficiency diminishes from 3% at 2^ = 1 GeV^ to 2% at = 10 GeV^. 
The data are corrected for this inefficiency. The systematic uncertainty on the trigger efficiency 
is estimated to be 1% for S9 and 0.5% for the other triggers. 

The inefficiency of the software filter farm component of the trigger is determined using a 
sample of the rejected events, recorded for cross checks. The primary reason for the rejection 
is the online reconstruction of the event vertex which occasionally wrongly classifies ep events 
as non-ep background. The loss corresponds to 0.7% for the NVX analysis and 0.8% (0.5%) 
for the SVX-BST (SVX-BDC) analysis. This loss is consistent with being uniform across the 
phase space and is applied as a global correction with a systematic uncertainty of size equal to 
the correction. 

The efficiencies of the electron identification requirements (cluster shape, hadronic fraction, 
BDC or BST validation) for high energies of the scattered electron are evaluated using events 
passing all other selection cuts but the one to be investigated. This direct approach is applicable 
for E'^ > 20 GeV due to the negligible background. 
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For low E'g, the electron identification efliciency is studied after the background subtraction. 
In addition, special background-free samples are used. One such sample comprises initial state 
radiation events with the radiated photon detected in the photon tagger. Background levels 
below 1% are achieved in this case, even for E[ ~ 3 GeV, by requiring the sum of the photon 
tagger and scattered electron energies to be close to the electron beam energy. The efficiency of 
the electron identification cuts does not vary much as a function of E'^. It is always above 90% 
and is well reproduced by the simulation. 

The efficiency to find a central tracker vertex for the BDC analysis is determined using 
events with a reconstructed EST vertex. As shown in figure 7, it is larger than 50% for y^. > 0.03 
and 93% for j/x > 0.1. The BDC analysis is restricted to i/s > 0.03. At larger values of y^. the 
efliciency decreases again, the effect being more pronounced in the data than in the simulation. 
The reason for this difference is a deficit of events with a large rapidity gap in the DJANGOH 
model, as already observed in [35]. This is accounted for by applying a correction to the MC 
simulation. The systematic error of the cross section resulting from this correction is found to 
decrease with increasing Q^, from 10% to 2%. 

A special procedure is developed to determine the EST inefficiency. Two main sources of 
inefficiency are distinguished which are both closely related to the readout procedure. The first 
is a hit finding inefficiency, which mostly depends on the performance of the front-end amplifier 
readout chip (APC). This efficiency is determined for each APC using EST tracks, requiring hits 
reconstructed in all but the EST plane under investigation. For most of the APCs the efficiency 
is high (> 95%), but about 5% of APCs have efficiencies below 80%. A few EST regions, with 
an APC efficiency below 40%, are excluded from the analysis. 

Correlated readout losses constitute the second source of inefliciency. In this case, signals 
are lost coherently in either ESTl, EST2 or in both modules. The main source of coherent 
losses comes from timing desynchronisation at a level of about 5% with some dependence on 
the (p sector. The coherent losses are measured separately and combined for ESTl and EST2 for 
each (p sector. They are measured using a background-free DIS sample (15 GeV< E'^ < 32 GeV) 
with a well reconstructed CT vertex and EDC track. All sources of EST losses are incorporated 
into the simulation. 

The efliciency of the EST track segment finder, used to reconstruct the event vertex within 
the nominal z range and to validate the SpaCal electron candidate, is checked globally for data 
and for the DJANGOH simulation using events with a well reconstructed central vertex and 
a EDC track segment. The photoproduction background is suppressed by demanding E[ > 
15 GeV. In this procedure not only is the EST efficiency examined, but the description of the 
EST acceptance and the imperfections of the tracking algorithm are also checked. Figure 8 and 
figure 9 show the global EST efficiency as a function of the electron candidate's radial position 
in the SpaCal for the NVX and SVX samples, respectively. Eased on this comparison, the 
systematic uncertainty attributed to the description of the EST efficiency is taken to be 2% in 
both the NVX and the SVX data analyses. This value also includes uncertainties arising from 
inefficiencies of the other electron identification criteria described above. 

The efficiency of the selection criteria based on the EDC, A^^iink bdc ^ 4 and |ArBDC-spaCail < 
1.5 cm, is determined for events with E'^ > 20 GeV for data and for the DJANGOH simulation. 
A correction to the simulated events is applied to account for discrepancies which are largest 
in the narrow transition region from small to large cells in the EDC. Events in this region are 
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rejected from the NVX-S9 analysis. The systematic uncertainty attributed to the BDC efficiency 
amounts to 1.5%, also covering differences between data and the MC simulation for other elec- 
tron selection criteria. 

8.2 Subdetector Alignment and Calibration 
8.2.1 Alignment 

Alignment of the Central Tracker, BDC and SpaCal. The relative alignment of the HI 
sub-detectors and the alignment of the detector with respect to the beam direction is performed 
in several steps. The first step is the internal alignment of the Central Tracker. The x and y 
coordinates are defined by the precisely known positions of the CJC wires while the z coordinate 
is defined by the COZ. Using cosmic muon tracks, the relative positions of the inner and the 
outer CJC parts, the location of the CIZ and the parameters for the z coordinate measurement in 
the CJC are determined. The beam axis is reconstructed by measuring the x and y coordinates 
of the interaction vertex as functions of its z coordinate. 

The alignment of the SpaCal and of the BDC is performed using high energy electron can- 
didates, with E'g > 20 GeV, linked to a central track. The central track is required to have at 
least two reconstructed CIZ hits and the 6 uncertainty must be smaller than 2mrad. The x 
and y offsets of the BDC and the SpaCal are measured by studying the difference in the polar 
angle measurements for these electron candidates between the Central Tracker and the BDC, 

= ^cT - ^BDc, and between the azimuthal angle measurements from the Central Tracker and 
the SpaCal, A^ = 0ct - ^spaCai, as functions of the azimuthal angle ^spaCai- The two methods 
find a consistent alignment in the x direction. For the y direction, the alignment is found to be 
different by 2 mm between the A0 and A0 methods. The average of the two values is used to 
correct for the misalignment. 

The z offset of the BDC is measured by studying A0 versus ^ct- The z offset of the SpaCal 
is checked by comparing the 6 measurements in the BDC and in the SpaCal. The tilts of the 
backward detectors are studied using A6' versus 6'ct for positive and negative x and y separately; 
they are found to be negligible. Figure 10 shows the comparison of the 6e measurement in the 
Central Tracker and the BDC after alignment. 

The SpaCal alignment with respect to the beam direction is cross-checked using quasi- 
elastic QED Compton (QEDC) events. These are scattering events of the type ep ^ epy with 
a hard photon radiated from the lepton line, the proton being scattered quasi-elastically at low 
momentum transfer such that the outgoing electron and photon are detected in the main detector, 
nearly back-to-back in azimuth. The QEDC process is selected by requiring two energy deposits 
in the electromagnetic SpaCal section with energies above 4 GeV. The sum of both cluster 
energies is required to exceed 25 GeV. The back-to-back requirement is enforced by demanding 
cos A0ey < -0.9 with A^g-j, being the azimuthal angle between the electron and the photon. 
Elastic events are selected by demanding no tracks reconstructed in the CJC and low activity 
in the calorimeters apart from the selected electron and photon. This alignment agrees within 
1 mm with the alignment obtained using central tracks. 

The dominant uncertainty of the alignment stems from the difference in the y direction 
between the A0 and A0 methods. Since the HI detector is nearly (p symmetric, biases in y 
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reconstruction do not lead to large shifts in the measured cross section. To cover a potential 
global bias of the Be measurement, a systematic uncertainty of 0.5 mrad is assigned for the polar 
angle reconstruction with the BDC and the Central Tracker. 

BST alignment. In the global BST alignment, the position of the BST is determined with 
respect to the HI coordinate system. In the internal BST alignment, radial offsets and rotations 
around the z axis of the individual wafers are determined. The global and internal alignments 
use the electron track reconstructed from the Central Tracker vertex and the BDC track segment 
as a reference and compare it to the track segment found in the BST. 

During the detector assembly each sensor is fixed to its nominal position with a mechan- 
ical precision of about 100/im. Remaining degrees of freedom are 128 radial shifts and 128 
rotations of the wafers. For the r strip sensors, these parameters are determined for all detec- 
tors simultaneously using the global minimisation package Millepede [84]. The degeneracy 
between shifts and rotations is resolved utilising the wafer overlap regions. Typical shifts are 
less than 200 jum and most rotations are less than 1 mrad. Figure 11a) shows the distribution of 
the number of BST linked hits as a function of (pe. Figure 1 lb) shows the difference in the Og 
measurement between the two BST overlapping sectors in these cases, after the BST alignment. 
An agreement to better than 0.2 mrad is observed. Based on this study, the uncertainty on the 
scattered angle reconstruction by the BST is taken to be 0.2 mrad. 

The alignment of the u strip detector is done in an analogous way. Here, shifts perpendicular 
to the u coordinate for the 8 wafers are determined simultaneously using the interaction vertex 
and the BDC measurement as additional external constraints. The shifts of up to about 100 yum 
are included in the external alignment. 

8.2.2 Electromagnetic Energy Calibration 

The largest uncertainty in the electromagnetic energy calibration stems from fluctuations of the 
gain factors of the individual SpaCal photomultiplier tubes. During the data taking, an initial 
cross calibration of the SpaCal cells was performed using cosmic muons. The stability of the 
gains was controlled by means of a dedicated LED system. First corrections to the gain factors 
were applied using DIS events based on the position of the "kinematic peak" - an enhancement 
in the E'^ distribution close to the electron beam energy which is characteristic of DIS at low 
at HERA. 

At the analysis stage, a cell-by-cell gain determination is performed using the double angle 
(DA) calibration. The DA method is also used to perform additional non- uniformity correc- 
tions taking into account variations of the energy scale on the sub-cell size level. The SpaCal 
energy non-linearity, caused particularly by the energy losses in dead material in front of the 
calorimeter, is modelled in detail [85] using the HI detector simulation based on the GEANT 
program [77]. The simulation is checked and corrected using ttP yy decays. Finally, the en- 
ergy scale is checked using Jlif/ ^ ee decays and QED Compton events. All calibration steps 
are described in the following. 
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Double Angle Calibration. The double angle calibration procedure makes use of kinematic 
peak events. Large statistics are available in this kinematic domain with negligible background 
contamination. For y < 0.1 the hadronic methods of y reconstruction (equations 13 and 14) have 
superior resolution. In particular, the scattered electron energy can be re-expressed in terms of 
the hadronic (equation 17) and electron scattering angles as 

EAl-ypA) tan^ 

^DA - g , ?/DA - ^ ^ • (^Uj 

sin^ -f tan ^ + tan -f 

In this method the scattered electron energy is calibrated to the electron beam energy^. The 
calibration corrects for genuine miscalibration and also energy loss in the dead material between 
the interaction point and the calorimeter. The same calibration procedure is applied separately 
to the data and the simulated events. 

For the calibration, events with E'^ > 20 GeV are selected. The event vertex position and 
the scattered electron angle are measured using the BST. A good ^da resolution is achieved by 
requiring 15° < dh < 80°. The calibration is performed by adjusting the gain factors of the 
individual SpaCal cells, such that the energy of the cluster agrees with the reference given by 
£'da- This is achieved in an iterative procedure: for each selected event, a ratio Cev = Eua/E'^ 
is calculated. The cluster energy is usually shared among several cells; the contribution of each 
cell with an energy Ei is given by V7, = EJE'^. A W,- weighted average of Cev for each calorimeter 
cell is then calculated based on all calibration events. This average is used to modify the gain 
factor at the next iteration. The calibration procedure converges after three iterations. 

The cell-by-cell calibration is followed by calibrations as a function of (i) the distance be- 
tween the centre-of-gravity of the cluster and the centre of the cell with highest energy to correct 
for biases of the clustering algorithm, (ii) 7?box = max(|jCsp|, \ysp\), where Xsp, y^p are the x and 
y cluster coordinates, in order to correct for energy losses in between SpaCal cells, and (iii) 

''sp = yj^lp + ?/sp' to correct for losses in the dead material in front of the calorimeter. These 
additional corrections are applied sequentially. 

The results of the double angle calibration are checked by comparing the distribution of the 
electron energy E'^ in the data and the simulation for the standard selection in the kinematic 
peak region. By comparing the widths of these distributions an additional Gaussian smearing of 
1.1% (0.2%) is applied to the electron energies in the simulated events for the NVX (SVX) data 
set. The need for this smearing in the MC may be due to short time scale drifts of the photo- 
multiplier gain factors which are not simulated, to imperfections in the shower shape simulation 
or to a deficiency in the passive material simulation. For the NVX sample, the kinematic peak 
comparison is presented in Figure 12a). Figure 12b) shows the Kolmogorov-Smirnov (KS) test 
probability distribution as a function of the relative shift in the energy distribution between the 
data and the simulation. Shifts above 20 MeV are excluded, which corresponds to a relative 
energy scale agreement better than 0.1%. The systematic uncertainty on the relative energy 
scale at the kinematic peak is taken to be 0.2% to account for the uncertainties of the HERA 
beam energy, for uncertainties in the resolution adjustment in the simulation, deficiencies of the 
double angle method and a residual variation of the level of agreement in the kinematic peak 
between data and MC for diff'erent bins. 

''The influence of initial state radiation, which eff'ectively reduces the electron beam energy, is small for this 
kinematic selection and is included in the simulation. 
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Calibration using n yy. The double angle energy calibration determines the energy scale 
of the individual SpaCal cells and radially dependent corrections of the energy loss for scattered 
electron energies close to the electron beam energy. The deviations from the linearity of the 
SpaCal response are measured using jip yy decays which probe much lower energies. 

Events with two clusters in the electromagnetic SpaCal section are selected. The larger of 
the two cluster energies is required to be above 2.2 GeV, exceeding the trigger energy threshold, 
the smaller cluster energy is required to be above 0.7 GeV. The event vertex is determined from 
tracks reconstructed in the Central Tracker. The two clusters are assumed to be produced by 
two photons. The invariant mass, Myy, is calculated using the reconstructed cluster energies and 
positions. 

The simulation of yy decays is checked using the PHOJET MC sample. A reasonable 
agreement of the simulation with the data is observed for the total energy of the two clusters as 
shown in figure 13a). The simulated ;r" energy spectrum is reweighted to that observed in the 
data, in order to reproduce the opening angle and individual photon energy distributions. 

The two-photon mass distribution is shown in figure 13b). A prominent peak above the 
background is observed close to the nominal tjP mass. The peak is shifted to lower values, 
around 130MeV instead of 135 MeV. This diff"erence is not reproduced by the MC simulation 
after the double angle calibration. The figure shows the distribution of simulated events after 
applying an additional correction of -3% to the energy scale for them. The data and the sim- 
ulation are then in a good agreement. The shift of the peak is possibly caused by not fully 
simulated energy losses in the dead material in front of the calorimeter. 

The low yy invariant mass and the relatively high photon energy cuts in the study of tjP yy 
decays lead to a rather small separation between the photon clusters in SpaCal, with an average 
separation of only 13 cm. An overlap of the adjacent clusters could lead to an energy scale 
shift. Additional studies are performed to estimate this eff"ect. The data sample is split into sub- 
samples with approximately equal statistics based on the larger or smaller cluster energy or on 
the cluster separation. In addition, the Myy distribution is studied as a function of the projected 
7r° location in the calorimeter, the latter being calculated as an energy weighted sum of the two 
cluster positions. In all these studies the relative shifts of the energy scale in the data versus the 
simulation are consistent within 1% which is taken as a systematic uncertainty of the energy 
determination at = 2 GeV. 

A check of the relative energy scale using decays is also performed for the SVX sample. 
The larger distance from the decay vertex to the calorimeter leads to larger separations between 
photon clusters, on average 18 cm. The relative shift of the Myy distribution between the data 
and the simulation after the double angle calibration is -2.7% in this case, consistent with the 
shift observed for the NVX sample. 

The relative bias of the energy scale is corrected in the data assuming a linear dependence on 
E'^. No correction is applied at E'^ = 27.6 GeV and a correction of -1-3% is applied dXE'^ = 2 GeV. 
The systematic uncertainty of the energy scale determination is also assumed to follow a linear 
dependence rising from 0.2% at E'^ = 27.6 GeV to 1% at 2 GeV. 

Tests of the SpaCal Energy Calibration. The SpaCal energy response is checked using 
Jjifj ee decays and QED Compton ep epy events. The y/i/f candidates are selected 
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by requiring exactly two electromagnetic clusters with a total energy of less than 22 GeV. At 
least one of the two clusters has to be linked to either a BST or a CJC track and both clusters 
must be associated with a BDC track segment. Events with additional CJC tracks not associated 
to the electrons are rejected, thus selecting events from elastic J I if/ production. The event vertex 
is defined by the CJC or the BST tracks. 

In this study, the SpaCal energy measurement is explicitly corrected to the absolute scale 
obtained from the mean ratio of the reconstructed to the generated electron energy from the 
DJANGOH simulation. Both the double angle and calibration corrections are applied, so 
that the peak in the di-electron invariant mass Mee distribution can be directly compared to the 
nominal Jjij/ mass, Mjj^ = 3.096 GeV. 

The distribution of for the NVX data is shown in figure 14. A clear enhancement around 
the nominal J/i/f mass is observed. The data are fitted with a sum of a Gaussian for the signal 
peak and a second order polynomial to describe the background shape. The fit uses the binned 
maximum likelihood method. The measured Gaussian peak position agrees with Mj/^ within 
1 .3 cr. Based on this agreement the deviation from the nominal energy scale is limited to be 
below 0.8% at 68% confidence level for energies of about 6 GeV. 

QED Compton events are used to check the calorimeter energy scale in the intermediate 
energy region. For elastic events, the energy of the scattered electron is related [86] to the polar 
angles of the scattered electron 6^ and the photon 6y by 



2Ee sin 6y 
sin 9e + sin 9y - sin (o^ + 



^QEDC ~ . . / \ • V^^-* 



The comparison of the measured electron energy with ^qedc '•^^'•^ '•^^ SpaCal energy scale 
linearity in the range 4-23 GeV. 

For the QED Compton energy scale check, a bias free reconstruction of the electron and 
photon angles is essential. Therefore in addition to the basic QEDC event selection described in 
section 8.2.1, both electron and photon SpaCal clusters are required to be linked to BDC track 
segments. This implies that the photon converted in the detector material in front of the BDC. 
The electron cluster is identified by requiring a BST link. The photon cluster must have no 
associated signals, neither in the BST nor in the CIP. 

The results of all calibration studies are summarised in figure 15. Both the J/ifr and the var- 
ious QEDC energy scale determinations are inside the uncertainty band. The scattered electron 
energy distributions and the uncertainty bands attached to the simulated E'^ distributions in the 
kinematic peak region are shown in figure 16 for the NVX and the SVX analyses. The data are 
well described by the simulations. 



8.2.3 Calibration of the Hadronic Energy Scale 

The calibration of the calorimeters employed for the hadronic final state energy measurement is 
based on kinematic constraints relating the scattered electron to the hadronic final state. For the 
calibration of the LAr calorimeter, conservation of the total transverse momentum Pj is used. 
The SpaCal calibration makes use of the conservation of E-P^. 
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Calibration of the LAr calorimeter. The hadronic final state in the central and forward re- 
gions of the HI detector is reconstructed using a combination of tracks and LAr energy deposits 
(see section 6.4). The LAr calibration coefficients are determined for the eight calorimeter 
wheels, each subdivided into eight octants, separately for the hadronic and electromagnetic 
sections. There are thus 120 calibration constants in total^, corresponding to the calorimeter 
segmentation in rapidity, azimuthal angle, and depth. The same calibration procedure is applied 
to the data and the MC simulation. 

To reduce the influence of the SpaCal on the calibration of the LAr, forward and central 
hadronic angles are selected: 13° < 9i, < 150°. The electron transverse momentum Pj is 
determined from the SpaCal energy and the 6^ measured by the BST. The photoproduction 
background is reduced to a negligible level by requiring E'^ > 20 GeV. 

In the calibration procedure, a least squares minimisation of the following function is per- 
formed 



(22) 



Here the transverse momenta Pj^ and P^ are the vector sums of the contributions from the 
tracks and the SpaCal, respectively, P;^ are the vector sums of the contributions from all cells of 
a calorimeter volume j, (ph is the azimuthal direction of the hadronic final state and aj are the 
calibration coefficients, which are free parameters. The outer summation is performed over all 
DIS events selected for the calibration. 

The Pj balance between the scattered electron and the calibrated hadronic final state is 
studied as a function of various variables, such as P^, 6^ and y^.- For central events, where i/^ > 
10"^, the simulation reproduces the behaviour of the data within 2% accuracy. At lowest y, the 
hadronic final state is produced at small polar angles and partially escapes the LAr acceptance. 
In this case, simulation and data agree within 10%. The systematic uncertainty of the hadronic 
scale is therefore extrapolated linearly in log y, from 10% at yj_ = 10"^ to 2% at j/x = 10"^. It is 
then set to 2% for j/s > 0.01. Figure 17 shows the overall Pt balance for the standard analysis 
selection. The vertical line at Pj/Pj = 0.3 indicates the analysis cut value. An increase in 
number of events for P'^/P^ < 0.3 corresponds to very low y < 0.001. The data agree with the 
simulation within the hadronic energy scale uncertainty. 



Hadronic Energy Calibration of the SpaCal. For large values of y > 0.4, the contribution 
of the SpaCal to the total E-Py becomes larger than the combined contribution of the LAr 
calorimeter and tracks. Given the accurate knowledge of the SpaCal linearity after calibration 
(section 8.2.2), the study of E-P, as a function of E'^ allows a check of both the linearity and 
the absolute scale of the SpaCal hadronic measurement to be made. 

The E-Py distribution is studied for E[> 1 GeV in E'^ intervals of 1 GeV. For each interval, 
a Kolmogorov-Smimov test is performed to estimate a possible shift in the E-Py distribution 
between data and simulation. For the SVX analysis, the data and the simulation agree well 
within their statistical uncertainties, while for the NVX data sample a global shift of ~ 1 GeV is 
observed. This shift is applied in the NVX analysis to the simulated events. An additional sys- 
tematic uncertainty, A(£'-i'^)spaCai = 0-5 GeV, is considered for both SVX and NVX analyses. 

^The most backward LAr wheel does not have a hadronic section. 
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Figure 18 shows the E-P^ distribution for the data and the simulation. The uncertainty band 
includes a +0.5 GeV variation of the SpaCal contribution to the total E-P^. 

8.2.4 Calorimeter Noise Uncertainty 

For y < 0.01, even a small fake energy contribution in the LAr can strongly affect the determi- 
nation of z//,. Therefore, a dedicated procedure is used to identify the LAr noise, as described 
in section 6.4. Samples of LAr electronic and beam induced noise are recorded in special runs 
and added to the simulation. 

The uncertainty of the noise influence on the DIS cross section measurement is determined 
as a function of yn by studying the ratio i//,,noise /(?//, + ?//i,noise), where ?//,,„oise is defined as z/A,noise = 
YiiiEi - P,j)/2Ee with the sum extending over the identified noise cells only. This comparison 
is shown in figure 19 for the NVX and SVX data samples together with contributions to yn 
from the tracks, LAr and SpaCal calorimeters. The noise fraction is described by the simulation 
within 10% accuracy which is taken as a systematic uncertainty. Note that at high y the noise 
fraction is small. More details on the LAr noise uncertainty estimation can be found in [82,83]. 

8.3 Background Subtraction 
8.3.1 Methods 

The dominant background source for this analysis arises from very low photoproduction 
events in which the scattered electron escapes detection in the backward beam pipe and a parti- 
cle from the hadronic final state mimics the electron. Other potential background sources arise 
from non-ep interactions. They are studied using non-colliding HERA bunches and are found 
to be negligible. 

For a fraction of photoproduction events the scattered electron is detected by the electron 
tagger of the luminosity system. These events are used to study the photoproduction back- 
ground. The acceptance of the electron tagger, which corresponds to the geometrical aperture 
of the detector as well as to the detection efficiency, is determined using Bethe-Heitler ep epy 
events [40], in which the scattered electron and the emitted photon are detected in the electron 
and photon tagger, respectively, and is parametrised as a function of y. The acceptance is large 
in the range 0.3 < y < 0.6. 

The simulated photoproduction background (PHOJET) is normalised based on events where 
the scattered electron is detected by the electron tagger and all of the analysis selection crite- 
ria^ are satisfied. Two normalisation methods are used. In the first method the background is 
normalised globally and then subtracted bin-by-bin 

Kp = Kmc- ^^'^' . <is=<ta-A^;.- (23) 



^For the electron method, this selection excludes the E-P- cut in order to increase the electron tagger accep- 
tance. In addition, to reduce the influence of overlapping DIS and Bethe-Heitler events, the absence of energy 
deposits in the photon tagger is required, and the total (fi-P )tot - E-P^ + lE^^^^^^, where ^fagger energy 
measured in the electron tagger, has to be less than 75 GeV. 
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Here, A^^j^ {N'^p) is the estimated number of DIS (photoproduction) events in the cross section 
measurement bin z, N^^^^^ and A^^^ are the numbers of data and PHOJET events in bin z, 
respectively, and A^tag, A^bg mc, tag are the total numbers of events detected using the electron 
tagger in the data and the PHOJET simulation, respectively. 

In the second method the background is normalised bin-wise using the bin-averaged tagger 
acceptance A, and then subtracted in each bin 

A^' A^' 

b" MC tas ; ■• tan 

= ' ^DIS = A^data - ^ . (24) 

where N^^^ and A^^^ are the numbers of events detected by the electron tagger in bin z in the 
data and the PHOJET event sample, respectively. Both methods lead to a cancellation of global 
selection uncertainties, while the second method (equation 24) also allows local uncertainties 
to cancel at the expense of an increased statistical uncertainty. 

For the NVX-S9 analysis, the global normalisation of the background (equation 23) is 
used, since for this sample the E'^ and 0^ distributions are well reproduced by the simulation 
(figure 20). Furthermore there is a direct control of the background normalisation as discussed 
in the next section. For the other analyses, a local bin-wise normalisation is performed. As 
a cross check, both normalisation methods are used for all samples, leading to cross section 
results consistent within statistical uncertainties. 



8.3.2 Normalisation Uncertainty 

The photoproduction background normalisation is checked for the NVX-S9 analysis using elec- 
tron candidates associated with tracks of opposite charge to the lepton beam charge, termed 
"wrong charge" tracks. Assuming charge symmetry of the background tracks, the wrong charge 
track sample gives an estimate of the remaining background in the correct charge sample. The 
track charge can be measured for tracks which are reconstructed in the BST sector equipped 
with u strip detectors in addition to the r detectors. 

In this method, any charge asymmetry creates a bias. In addition, the requirement of a w 
strip track in the background study could modify the normalisation compared to the standard 
sample. The geometrical acceptance and efficiency e of the u strip track reconstruction are first 
determined based on a high E'^ sample in which the background can be neglected. The accep- 
tance and efficiency are well described by the simulation. The acceptance difference between 
data and simulation is found to be (2.0 + 1.3)%, while the efficiency difference is determined to 
be (0.2 + 0.5)%. 

All events within the u sector acceptance passing the NVX-S9 analysis cuts, A^acc, are clas- 
sified according to A^acc = Nq + N-i- + N-, where A^o denotes all events without a linked u track, 
A'^+ is the number of all events with correct sign tracks (positive, as expected from the scattered 
positron) and A^_ is the number of all events with wrong sign tracks. If k = N^^/N^^ is the 
background charge asymmetry ratio, then the total number of background events in the u sector 
geometrical acceptance is 

= A^o + A^-(l + ^) + A^+|l - 



29 



Here A^acc(A^sig) denotes the number of accepted (genuine DIS signal) events. 

The charge asymmetry of the background for the PHOJET simulation is found to be /c = 
0.79+0.06. A dedicated study of the origin of this asymmetry [37] showed that the main effect is 
due to the difference between the proton and antiproton interaction cross sections and the visible 
energy which they deposit in the SpaCal. A larger value of \Elp\ is expected for antiprotons 
since they annihilate at the end of their paths. Indeed, for simulated events with \Elp\ > 2 the 
deviation of k from unity is larger: k = 0.60 ± 0.14. From the data with \E/p\ > 2 a consistent 
value K = 0.65 + 0.12 is measured*^. The charge asymmetry is also checked using events in 
which the scattered electron is detected in the electron tagger. It is found to be 0.82 + 0.17. 
The PHOJET based asymmetry estimate is also consistent with the value estimated in [37] 
using tagged events, k = 0.91 + 0.04. In order to cover the findings on the charge asymmetry 
explained above, a value k = 0.9 + 0.1 is assumed for this analysis. 

The ratio of the number of photoproduction events obtained using equation 25 to the es- 
timated number of events based on the electron tagger, equation 23, for the E'^ range of the 
NVX-S9 analysis, is r = 1.00 ± 0.14stat ± 0.05asym- Here the first error gives the statistical 
uncertainty and the second error corresponds to the uncertainty in the background asymmetry 
determination. Figure 21 shows the distribution of E'^ for the background events, estimated 
using u sector tracks. The systematic uncertainty on the background normalisation is taken to 
be +15%, based on the statistical uncertainty of the u sector sample and the uncertainty in the 
background charge asymmetry. 

8.4 Luminosity Determination 

The luminosity measurement is based on Bethe-Heitler events detected using the photon detec- 
tor. A precise luminosity measurement requires a good understanding of the beam optics, of the 
photon detector acceptance and its variation with changing beam conditions. The uncertainties 
related to the acceptance are similar for the NVX and the SVX data. 

The time structure of the ep bunch crossings is characterised by the main proton bunch 
accompanied by satellite bunches. Two such bunches are at +4.8 ns away from the nominal 
bunch and lead to ep interactions at about +70 cm from the mean vertex position. The photon 
detector is sensitive in a time window of about 12 ns for Bethe-Heitler events and thus does 
not distinguish interactions at the nominal vertex position from satellite bunch interactions. 
The luminosity measurement therefore requires the fraction of satellite bunch interactions to be 
determined independently. This is possible in HI using TOE and PPU systems. 

Eor the SVX data, with the main bunch centred at z = 70 cm, the backward satellite is 
located at z ~ cm. The backward satellite in this case gives a larger contribution to the lu- 
minosity measurement than the forward satellite at z ~ 140 cm. The fraction of events in the 
backward satellite can be determined directly from the fraction of DIS events with a recon- 
structed vertex around z = cm and amounts to 2.7%. A 3% uncertainty is assigned to the 
luminosity measurement for the SVX data, which covers the differences observed between the 
methods of determining the satellite bunch fraction and also includes uncertainties related to the 
photon detector. The same procedure is performed to verify the contribution from the forward 

**At low energy, the contribution of DIS electrons with \E/p\ > 2 is negligible. 
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Uncertainty 



BST efficiency 2% (BST) 

BDC efficiency 1 .5% (BDC-Central vertex) 

Central Tracker vertex efficiency 2 - 10% (BDC-Central vertex) 

Trigger efficiency 0.9% (NVX) 



1.1% (NVX-S9) 
0.9% (SVX-BST) 
0.7% (SVX-BDC) 
Radiative corrections 0.5% 



Table 4: Summary of the systematic uncertainties. For the correlated sources, the uncertainties 
are given in terms of the uncertainty of the corresponding source. The effect on the cross section 
measurement varies from bin to bin and is given in table 10-14. For the uncorrected sources, 
the uncertainties are quoted in terms of the effect on the measured cross section directly and the 
type of analysis is given in brackets. 

satellite at -1-70 cm of the NVX data sample. In this case the different methods are in agreement 
within 0.7%, leading to a total luminosity uncertainty of 1.1%. 

In the course of this analysis an extended reanalysis of the 1997 data at < 12GeV^, 
this sample termed B in [37], was performed, which reproduced the published cross sections in 
shape. These, however, are to be multiplied by a factor of 1 .034 as the result from an improved 
analysis of the satellite bunch structure and the photon detector acceptance. This corresponds 
to a shift of two standard deviations of the quoted luminosity measurement accuracy. 

8.5 Summary of Systematic Uncertainties 

The systematic uncertainties are classified into two groups, bin-to-bin correlated and uncorre- 
cted systematic errors. For this analysis the correlated sources are the electromagnetic and 
hadronic energy scales, the electron scattering angle, the calorimeter noise, the background 
subtraction and the normalisation uncertainty. The uncorrelated errors are related to various 
efficiencies and radiative corrections. A summary of the correlated and uncorrelated errors for 
the present analysis is given in table 4. 
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The large overall contributions to the total error are due to the BST electron track recon- 
struction efficiency and the Central Tracker vertex efficiency uncertainty. The correlated error 
sources affect the DIS cross section measurement in a manner which depends on the kine- 
matic domain. The most pronounced variation arises with the inelasticity y. For high y, the 
uncertainty is dominated by the photoproduction background (about 6% for y = 0.8). For in- 
termediate y ~ 0.1, the E'^ scale uncertainty becomes more prominent for the electron method 
(about 3% cross section uncertainty). For y < 0.01, the dominant error source is the LAr noise 
(up to 10% cross section error). 

8.6 Control Distributions 

Data and Monte Carlo simulation distributions of important quantities for the events passing 
all selection criteria are compared in figures 22 - 25. Only events corresponding to analysis 
bins passing the stability and purity criteria are considered. The simulated distributions are 
normalised to the measured luminosity. The DIS MC cross section prediction is reweighted to 
a parameterisation using the fractal model introduced in section 11.1. A rather good (NVX) to 
acceptable (SVX) overall agreement is obtained in the description of the data by the simulation. 

Figure 22a-d) shows basic kinematic and vertex distributions for the NVX analysis. The 
background from photoproduction events is very small. It is larger at lower scattered positron 
energies E'^ as can be seen in figures 22 e) and f), which show the E'^ and 6^ distributions for the 
dedicated high y analysis (NVX-S9). In figure 23 basic kinematic distributions for the SVX- 
BST analysis a)-c), the SVX-BDC analysis d) and the SVX-BST analysis considering events 
from ISR bins only e), f) are shown. The ISR distributions are very well reproduced by the 
simulation. The other SVX plots reveal a small normalisation difference. Figure 24 shows the x 
and distributions for the two kinematic reconstruction methods, electron and S, in the NVX 
analysis. Figure 25 shows similar distributions for the SVX analysis. Events are only taken 
into account from bins which pass the stability and purity criteria and are covered by the chosen 
method. For the SVX sample the data are less well described than for the NVX sample, but 
consistency is observed within the total measurement uncertainty including a 3% normalisation 
error of the SVX data. 



8.7 Cross Checks 

The stability of the cross section measurement is tested with a set of dedicated cross checks 
which can be divided into three classes: (i) checks for a given data set and a given reconstruction 
method, (ii) checks of the consistency between the different reconstruction methods, and (iii) 
checks of the consistency between the different data sets. 

The consistency of the cross section measurement for a given data set (e.g. NVX) and a given 
reconstruction method (e.g. the electron method with 6^ measured by the BST) is studied by 
splitting the data into two approximately equal sub-samples and comparing these sub-samples 
to each other. The data are compared as measured with the upper and the lower half of the 
SpaCal, for negative and positive z- vertex positions, and dividing the sample into an early and 
late data taking period. These tests are sensitive to local effects like efficiency variation, energy 
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miscalibrations and the stability of the luminosity measurement. In such studies no significant 
deficiencies in the data are observed. 

For the comparison of the cross section measurements for a given data set but using differ- 
ent reconstruction methods, the test samples are strongly correlated. The uncorrected statistical 
uncertainty is estimated in this case by subdividing the simulated events into a number of in- 
dependent sub-samples of equal size. The measurement of the cross sections is repeated for 
each sub-sample and the statistical uncertainty is calculated as the luminosity rescaled RMS of 
the resulting variations of the cross section measurements. Employing this technique, the cross 
section measurements based on different triggers and different 9e reconstruction methods (BDC 
and BST) are compared. In most of the cases the measurements with each of the samples agree 
within the uncorrelated statistical uncertainty. In a few cases the measurements agree within the 
total uncertainty only. A particularly interesting test is the comparison of the cross section mea- 
surement performed with the electron and sigma methods, since the two methods have different 
sensitivities to systematic error sources. The two methods can both be applied in many common 
bins where the purity and stability of the measurement are high for both methods. Figure 26 
shows an example of this comparison, performed for the NVX-BST data set. 

The third class of cross checks compares the cross section measurements performed with 
different data samples: SVX is compared to NVX and the new data are compared to the pre- 
viously published results. This comparison is an integral part of the cross section averaging 
procedure, as discussed subsequently. 

8.8 Cross Sections 

The cross section data measured from the SVX and NVX data samples are given in tables 10-14 
and presented in figure 27. The uncertainty of the new data is typically 3 - 4% and larger at 
the acceptance edges. Lowest values of Q^, down to 0.2 GeV^, are reached with the shifted 
vertex data. The analysis of the SVX data is mainly based on the BST but complemented by 
an independent analysis using the BDC at lower radii. For between 0.5 and 3.5 GeV^, the 
NVX and SVX cross section data overlap in their kinematic coverage and are observed to be in 
agreement. The kinematic region of larger > 5 GeV^ is covered by the nominal vertex data. 
The data at highest y, corresponding to smallest x, are obtained using the dedicated trigger S9 
and can be seen in figure 27 to be consistent with the behaviour of cr,. towards small x. 



9 Combination of HI Cross Section Measurements 

The new data cover a kinematic region which overlaps with data sets taken at 820 GeV proton 
beam energy in 1995 [35] and in 1997 (sample B) [37]. The combination of all these data, 
as described subsequently, provides a single data set in the range 0.2 < < 12GeV^ and 
5 • 10"^ <x< 0.02. 
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9.1 Procedure 



The combination of the data sets is based on the prescription introduced in [87] which is ap- 
plicable if the uncertainties of the measurements do not depend on the central values. This 
procedure is described in section 9.1.1. For the cross section measurements the estimated sta- 
tistical and systematic uncertainties depend on the central values. This leads to a modification 
of the averaging procedure as is described in section 9.1.2. 



9.1.1 Linear Averaging 

The averaging procedure is based on a x'' minimisation. For a single data set, the function 
can be defined as 



i I j «j 

Here /u' is the measured central value at a point / with combined statistical and uncorrelated sys- 
tematic uncertainty A,- = {A^^^^i + ^funcor)''^^- Further, aj denotes the central value determined for 
a correlated systematic error source of type j with an uncertainty Aq,., while d/u' /daj quantifies 
the sensitivity of the measurement /u' at the point / to the systematic source j. The function ;^fg^p 
depends on the set of underlying physical quantities m' (denoted as the vector m) and the set of 
systematic uncertainties aj (a). For the reduced cross section measurements one has fj.' = a\., 
i denotes a (jc, Q^) interval, and the summation over j extends over all correlated systematic 
sources. 

Introducing the variables bj = (aj - aj)/Aaj and F^ = (dju'/daj)Aaj, equation 26 can be 
written as 

i ' j 

If several analyses provide a number of measurements at the same (x, Q^) values, they can be 
combined using the formula above, generalised for the case of multiple data sets. Then a total 
function, is built from the sum of thexlxp functions for each data set according to 

Nm Im'-Y'^s p U._i.if Ns 

EE ^ ^'y-' '^'^ ^..-^Ef-;. (28) 

e i=l i,e j=l 

where the summation over i (j) runs over all Nm measured points (all A^^ systematic error 
sources) of all data sets considered. The symbol wie is equal to one if data set e contributes a 
measurement at the point i, otherwise it is zero. Similarly, the symbol Y'j^ equals to zero if the 
measurement i from the data set e is insensitive to the systematic source j. This definition of 
Xtot assumes that the data sets e are statistically uncorrelated. The systematic error sources bj, 
however, may be either uncorrelated (separate sources) or correlated across data sets (diff"erent 
data sets sharing a common source). 



34 



Since Xtot ^ quadratic form of m and b, it may be rearranged such that it takes a form 
similar to equation 26 



1=1 

Ns Ns 



Ah 



(29) 



y=l A=l 



The data averaging procedure, described in detail in appendix A, determines the average values 
/u''^^^, the uncorrected uncertainties A,;ave, the average systematic error source values ay.ave, the 
dependencies of //''""'^ on aj, d/j''""" /daj, and the matrix (A'g)jk. The value ofx^^j^ corresponds 
to the minimum of equation 28. The ratio x^i^/ '^dot is a measure of the consistency of the 
data sets. The number of degrees of freedom, n^^ot, is calculated as the difference between the 
total number of measurements and the number of the measured points Nm- This procedure 
represents a method to average data sets, which allows correlations among the measurements 
due to systematic uncertainties to be taken into account. 

The matrix (A^ )jk can be diagonalised and the;^;^ function takes a form similar to equation 27 



Ns 



1 ^ j "j 



i,ave 



+ 



Ns 



(30) 



!=1 



7=1 



where b'j = 2a- Ujk{bk - fik,ave)Djj and fik^ave = Q^/t.ave/Ao.^.. The orthogonal matrix U connect- 
ing the systematic sources before and after averaging and the diagonal matrix D are given in 
appendix A. 



9.1.2 Implementation for the Cross Section Averaging 

The function of equation 26 is suitable for measurements in which the uncertainties are 
absolute, i.e. do not depend on the central value of the measurement. However, for the HI cross 
section data considered here, the correlated and uncorrected systematic errors are to a good 
approximation proportional to the central values (multiplicative errors), whereas the statistical 
errors scale with the square roots of the expected number of events. In this case the combination 
of the data sets using equation 26 leads to a small bias to lower cross section values since the 
measurements with lower central values have smaller absolute uncertainties. To take this effect 
into account, the;^^^ definition is modified to 

' ^tsm^' ['^' - ^7 yj'^'^j) + uncor m'j j 

Here Yj = ^'j/fJ^', ^r,stat = A,-_stat//^' and 5,-,uncor = A/,uncor//^' are relative correlated systematic, 
statistical and uncorrected systematic uncertainty, respectively. This x^ definition is used for 
the averaging procedure and also for the phenomenological analysis of the data (see section 10). 
Equation 31 is equivalent to the one used in previous HI analyses [37] up to modifications in the 
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denominator. In contrast to equation 27, thex^ function of equation 31 is not a simple quadratic 
form with respect to m' and bj. The average is found in an iterative procedure: first equation 27 
is used to get an initial approximation for fi''^"^ and ySy.ave which are used to recalculate the errors 
asr^. = r^yu'-^™ and = 5^,t,tju'(//'-"^'=-2jr>''''™>S^-,ave) + (5,-,uncor/i''''''^ Then the determination 
of fj.^'^"^ is repeated. Convergence is observed after two iterations. 

For measurements with multiplicative errors the geometric mean instead of the arithmetic 
mean can be used as an alternative, i.e. the average is performed for Incr).. In this case the 
quadratic equation 27 can be used by replacing //' In cr[. A, (<^^stat "'"^funcor)^^^ ~^ y)- 

This logarithmic averaging procedure is used as a cross check. 

For the NVX and SVX analyses, the measured cross section values cr'., the statistical and 
uncorrected uncertainties 5,,stat, 5,,uncor and all correlated systematic uncertainties y'^ as well as 
the total error 5^oi = [<^fstat + ^^uncor + Z/t})^]'^^ are given in tables 10-14. The average of the 
HI data is reported in tables 15-20, where the average reduced cross sections cr^f^'^ = fi''^^^, the 

statistical (5,-,ave,stat, uncorrected 5,-,ave,uncor, correlated yf^'' and total 5,-,ave,tot = [^^ave,stat+<^^ave,uncor + 

J^jiyY''^)^V^^ uncertainties are given. The transformation matrix U is given in table 21. The 
shifts of the central values of the systematic error sources, in units of the systematic errors 
aj^ave/^aj, ^rc givcn in table 5. 

9.2 Compatibility of SVX and NVX Data 

The combination of the SVX and NVX data depends upon assumptions on the correlation be- 
tween different data points, within a data set as well as across the data sets. For each data set, 
two types of systematic uncertainty are considered: fully correlated ones, which are treated as 
aj in equation 31, and fully uncorrelated ones, which are added to the statistical uncertainties in 
quadrature and treated as Sj in equation 3 1 . Following the notation in table 4, the six sources of 
correlated uncertainties are E'^ scale, 6^, LAr hadronic energy scale, LAr noise, SpaCal hadronic 
scale and photoproduction background. A further correlated uncertainty arises from the lumi- 
nosity measurements. Concerning the relation between data sets, the systematic uncertainty of 
the luminosity measurement is separated into a 0.5% fully correlated theoretical uncertainty and 
an uncorrelated experimental part due to variations of beam and detector acceptance conditions. 
The other systematic uncertainties are considered to be uncorrelated. 

The systematic uncertainties which are correlated between data points can be assumed as 
either fully correlated, uncorrelated or partially correlated between the NVX and the SVX data. 
The reasons for correlations between data sets are the similarity in the calibration procedure and 
the detector setup. Uncorrelated effects arise from variations with time, differences between the 
kinematic ranges of the calibration samples, the dead material, the detector illumination or the 
acceptance. For each source the uncorrelated part is more important and thus all sources are 
considered to be uncorrelated between the NVX and the SVX data. 

To check the sensitivity of the averaged result to the correlation assumptions, the average of 
the NVX and SVX data, obtained by considering the six systematic sources to be uncorrelated, 
is compared to 2^ - 1 other possible assumptions in which each source is either fully correlated 
or fully uncorrelated. Most of the cases lead to numerically small variations for both the central 
values and the total errors of the average data. The only significant variation is observed for 
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the lowest y points for > 2 GeV^, if the LAr noise is assumed to be correlated between the 
NVX and SVX data. Since the LAr noise, however, is a time dependent uncorrected source, 
no additional systematic uncertainty is assigned to the combined measurement. 

The NVX and SVX data sets are fully consistent, according to the averaging procedure, 
with;^'^!^/ ^dof = 19.5/39. The shifts of the central values of the systematic uncertainties do not 
exceed one standard deviation. 

9.3 Global Combination of Low HI Data 

The new data given in tables 10-14 are combined with the previously published HI data obtained 
for a similar kinematic region. The comparison of the present cross section data, obtained 
by averaging the SVX and the NVX data, with the published cross section data, is given in 
figure 28. The new data are in agreement with the published NVX97 data [37] taking the 
+3.4% normalisation shift of the published data (section 8.4) into account. The data are also 
consistent with the SVX95 data [35] within their rather large uncertainties. For the combination 
of all data, the systematic uncertainties are considered to be uncorrelated across the data sets. 

The published HI data [35,37] were taken with a proton beam energy of Ep = 820 GeV 
Therefore a centre-of-mass energy (CME) correction, based on equation 1, is applied when 
comparing to the previously published cross section according to 

cr',^\x, Q^) = af\x, (£-) + Q") [/(i/82o) - /(?/92o)] ■ (32) 

Here crp"(x, 2^) is the reduced cross section rescaled to Ep = 920 GeV; crf'^{x, Q^) is the 
measured reduced cross section for Ep = 820 GeV; i/82o and 2/920 are the inelasticities for the 
two proton beam energies calculated as y = jAEeEpX, and Q^) is calculated using the 
fractal model for F2{x, Q^) and R = 0.5. This correction becomes large only at high y. To 
avoid any sizeable effect of the energy dependence of cr,- on the combination of the 820 and 
920 GeV data, the combination of the points at the same (x, Q^) is restricted to a region of 
?/82o < 0.35. At higher y the measurements are kept separately but they are affected by the 
combination procedure. The residual dependence on the Ft assumption for the average points 
is negligible. For illustrative purposes, the CME correction is applied to all 820 GeV data points 
in figures 28-32. 

The HI data sets are consistent with each other. If all samples are averaged in a single step 
one obtains ;^^ot/ ^dof = 86.2/125. Shifts of the central values of the systematic sources aj/Aaj 
are given in table 5. The systematic shifts imposed by the averaging procedure are mostly within 
one standard deviation. The most noticeable effects are a downward shift of the normalisation 
of the SVX95 data and a modification of the LAr hadronic energy scale of the SVX data which 
corresponds to a small adjustment of the SVX data at large x. 

The combination of the HI data using the;^'^ definition of equation 3 1 has been compared to 
that using the;^'^ definition of equation 27 and also using the logarithmic averaging procedure. 
For the bulk of the phase space, the definition of equation 27 would lead to a change of typically 
-0.7%, which increases to -2.0% for the data at < 0.5 GeV^. For the logarithmic average 
the difference compared to equation 31 is typically below 0.1%. 
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Table 5: Shifts of the central values of the systematic uncertainties aj^ave/Saj based on the aver- 
age of the published Ep = 820 GeV and the new NVX-SVX data. For example, the quoted value 
for the luminosity shift of the SVX95 sample, -1.60, corresponds to a -1.60 x 3% = 4.80% 
downward shift of the SVX95 cross section values. 

9.4 Combined Cross Section Results 

The combined low data and the resulting uncertainties are listed in tables 15-20 and shown 
in figures 28-32. There are 149 data points. The total uncertainty in the central Q^,x region 
of this measurement is about 2% but it becomes larger towards the edges of the covered phase 
space. At high y, for example, the measurement at a value of a few GeV^ has an uncertainty 
of about 5%. 

Figures 29 and 30 show the combined HI reduced ep cross section measurement and dif- 
ferent phenomenological descriptions as introduced below. For all bins, starting at large x 
the reduced cross section first increases for x ^ 0. For > 0.6 GeV^ there is a characteristic 
turn over of the cross section observed at the smallest x values. This region, for each Q^, corre- 
sponds to highest inelasticity, y = Q^/{sx), and thus the turn over at y x 0.6 can be attributed to 
the influence of the longitudinal structure function Fi. 

For y < 0.6 the influence of the longitudinal structure function is small and thus one can 
extract the structure function F2 with only a small residual dependence on the assumption on 
Fl. Using R = 0.5, F2 is extracted and shown in figure 31. The structure function F2 exhibits a 
steady increase as jc ^ for all bins. 

Figure 32 shows the measurement of the virtual photon-proton eff'ective cross section cr'^^^ 
as a function of at various values of W. The HI data are compared to the data of ZEUS [38, 
39] and to different models, as discussed below. A good agreement between the data sets is 
observed. The HI data extend the HERA measurements to higher and lower W and also cover 
the ~ 1 GeV^ region. 

10 Cross Section Analysis 

10.1 Rise of F2 at Low x and Extraction of R 

The rise of the structure function F2 towards low x has previously been described by a power 
law in X, F2 = c{Q^)x''^^^~\ where the exponent A increases approximately logarithmically with 
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In for 2^ > 2 GeV^ [88]. This simple parametrisation has been shown to model the ep data 
well for ;c < 0.01. 



This idea can be extended to fit the reduced cross section cr,. in order simultaneously to 
extract the exponent X and to estimate the longitudinal structure function F^. The measured 
ep cross section is sensitive to the longitudinal structure function only for large y > 0.5, a 
region which corresponds to a limited x range for a given value. Gluon dominance at low 
X suggests that the function may exhibit an x dependence similar to F2. In the subsequent 
studies using this ansatz it is assumed that Fl is proportional to F2 and that the coefficient of 
proportionality depends only on Q^. For the extraction of Ft, the ratio of ctl/ctt = R is used 
such that 

Fl(Q\ X) = F2(Q\ X) Y^J^y (33) 

and 

r p/'n2 

crr{Q\x) = c{Q^)x-'^Q^ 



1 - /(?/)- 



1 +R{Q^)\ 



(34) 



The combined 1995-2000 HI low data are fitted following equation 34 for each bin. 
These fits describe the data very well, as is illustrated in figure 29. The results of these fits 
are shown in figures 33 and 34. The fit results for A{Q^) are given in table 6. One can see in 
figure 33b) that the parameter A shows an approximately linear increase as a function of In 
for > 2 GeV^ as has been observed previously [88]. For lower the variation of A is 
diminished but relatively large uncertainties prevent definite conclusions. The normalisation 
coefficient c{Q^) rises with increasing for < 2 GeV^ and is consistent with a constant 
behaviour in the DIS region, as in [88]. 

The values of the coefficient R(Q^) are consistent with no dependence on Q^. The mean R 
is 0.55 + 0.05 with^ x^/ ^dof = 7.9/(8 - 1). While the experimental error is small there is a very 
strong model dependence, different parameterisations for F2 leading to significant changes in 
Fl, see section 1 1 . The value of the average R obtained in this model is consistent within about 
one standard deviation with R = 0.5 or ctl = \(Tt- This value of R leads to an Fi which is 
higher than the first direct measurement of Fi at low x performed by the HI collaboration [89]. 
The data in [89] correspond however to higher values (> 12 GeV^). 



10.2 Determination of using the y Dependence of the Cross Section 

The turn-over of the measured DIS cross section for the highest y values, apparent in figure 29, 
can be used for an extraction of the longitudinal structure function Fi using the so-called deriva- 
tive method [37]. The derivative of the reduced cross section with respect to In y is 



dcr,- 



dln i/ 



dF2 2y\2-y) ^ y^ dF^ 
- . . ...... • (35) 



dlnx (l+(l-?/)2)2 \+{l-yfd\nx 



At high y for a wide variety of models the term proportional to F^ becomes numerically larger 
than other contributions. Therefore the extraction of the derivative provides means for deter- 
mining Fl at low x and with little phenomenological assumptions. 



^For the determination of the mean, R values from different bins are assumed to be uncorrected and total 
errors are used. 
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0.003 


0.005 
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0.007 


8.50 


0.216 


0.003 


0.005 


0.003 


0.007 



Table 6: Results of the fit (equation 34) to the combined HI low data on the exponent A 
with the statistical ^.j^stat, uncorrelated systematic (5i,uncor, correlated systematic ^.i^cor, and total 
uncertainties 5^_tot- 

Experimentally, dcr,7d In y is approximated by i/avAcr,./Ai/, which is calculated for each pair 
of cross section measurements in neighbouring bins. Here Ay is determined using bin centre 
values, and y^y is the logarithmic average value. Only the Ep = 920 GeV data are used in this 
determination. The HI data are illustrated in figure 35 and are compared to the fractal and dipole 
models discussed below in sections 11.1 and 11.2. Similar analysis using the Ep = 820 GeV 
data was presented in [37]. The systematic uncertainties are evaluated by changing the cross 
sections for each source of systematic uncertainty and repeating the calculation of the cross 
section difference. For the model predictions, AcrJAy is calculated in an analogous way and 
using the same binning as for the data. 

For the extraction of the structure function Fl, the fractal fit, introduced in section 11.1, is 
used to estimate the dF2/dlnx contribution to Act,- I Ay, and also for the bin centre correction. 
To reduce the dependence on F2, only AcTrlAy value corresponding io y = 0.735 are used to 
determine Fi. 

The resulting longitudinal structure function values are shown in figure 36. The derivative 
method is only weakly dependent on the model assumptions. There are however large exper- 
imental uncertainties, mostly due to statistics and the photoproduction background at large y. 
The Ft data are consistent with a constant R = 0.5, as introduced above, and also with smaller 
values on R, as obtained in the dipole models. The dependence of the measurement on the 
assumption made for F2 is estimated by a comparison with results obtained when assuming F2 
to be independent of x. The difference between the extracted Ft values is shown as the band at 
the bottom of figure 36. 



11 Model Comparisons 

In the following the combined data are analysed in the context of the fractal model [33] and two 
versions of the colour dipole model [42,45], which unlike pQCD may be applied to describe the 
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Ql (GeV^) 


0.093 


0.010 




0.56 


0.07 



Table 7: Parameters of the fractal fit and their total uncertainties. For the central fit D2 is kept 
constant: D2 = 1.08. If the parameter Dj is floated, the fit gives D2 = 1.061 + 0.012. 

transition region from photoproduction to deep inelastic scattering. Fits are performed using 
equation 31. 



11.1 Fractal Fit 



In the fractal ansatz [33] , the proton structure function F2 is parameterised using five parameters 
Qo and Do to Di, as 

X i + ^ -1 ■ (36) 

The parameters of this model are determined with a fit to the cross section data, except for the 
parameter D2, which governs the structure function behaviour for the photoproduction regime 
and is fixed to D2 = 1.08. This parameterisation is used in the Monte Carlo reweighting pro- 
cedure. The fractal model [33] does not provide predictions for Fl. The same prescription is 
followed as for the A parameterisation fit described in section 10. 1 taking the Fl contribution to 
be proportional to F2. 

The values of R are found to be consistent with the A fit and with being independent of 
Q^. Thus for the fractal parameterisation of the reduced cross section, R is taken to be a 
constant, which results in the simple five parameter representation used in the present anal- 
ysis. The parameters of the fit are given in table 7. The fit describes the data well with 
X^/ ^dof = 155.3/(149 - 5). Similarly to the A fit, the value of R = 0.56 + 0.07 is consis- 
tent within about one standard deviation with R = 0.5. This agreement with the A fit may be 
attributed to the structure function F2 having a power law-like x dependence. 



F2iQ\x) = DoQl 



-D2+1 



I+D3- Di \nx 



11.2 Dipole Model Fits 

In the GBW model [42] the dipole-proton cross section & (see section 3) is given by 

&(x, r) = (TO |l - exp [-rV (^rl(x))]} , (37) 

where r corresponds to the transverse separation between the quark and the antiquark, and is 
an X dependent scale parameter, assumed to have the form 

rl(x) - (x/xoY . (38) 
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UIlCvl 1*11111 J' 


CT (mb) 


OA ^ 


U.J 


A 


0.256 


0.003 


Xq 


0.60x10-4 


0.03x10-4 



Table 8: Parameters of the GBW dipole fit and their total uncertainties. 



Parameter 


Value 


Uncertainty 


Rum (fm) 


0.605 


0.008 


A 


0.260 


0.003 


Xq 


0.45x10-4 


0.03x10-4 



Table 9: Parameters of the IIM dipole fit with A^o = 0.7 and their total uncertainties. 

For small r tq, & is proportional to (colour transparency, & ~ {rllrof) while for r » 
tq the cross section approaches a constant value (saturation, & ^ ctq). The boundary in the 
(x, Q^) plane which separates these regions is described by the "critical line" at the x dependent 
saturation scale Q^^^x) = \lr^{x). The model provides predictions for both ctt and ctl in terms 
of only three parameters, ctq, Xq and A. 

The fit to the reduced cross section with the dipole model of GBW ("GBW fit") yields a 
I ^doi = 183.1/(149 - 3), acceptable but worse than that for the fractal model. It has been 
suggested that improved models of & lead to a better description of the data and a variety of 
models has been developed. As an example, a fit using & as proposed in the IIM model, with 
A'^o = 0.7 as defined in [45], has been performed. This fit also has three free parameters and gives 
I f^dof = 178.2/(149 - 3). The results of the two dipole model fits are shown in figures 30-32 
and 35-36. The dipole model fit parameters are given in tables 8 and 9. 

To trace the origin of the differences between the models, predictions for the structure 
functions F2 and Fi are compared individually. As an example, figure 37 shows the comparison 
between the three models for the bin = 1.2 GeV^. The structure functions F2 agree rather 
well for the models considered for x > = 0.18 x 10 ^, where x^ corresponds to the saturation 
radius of the GBW dipole model at the chosen value. However, for x < x^ the dipole models 
show a softer F2 dependence on x. This holds in particular for the IIM dipole model. The 
main difference between the models is in the structure function F^. As shown in figure 37, the 
predictions of the dipole models are nearly half of the result for F^ obtained with the fractal 
model analysis. 

The strict correlation between F^ and F2 predicted by the dipole models could be broken 
by higher twist effects [90]. To quantify the influence of the structure function Fi another fit to 
the reduced cross section data is performed, in which the F^ prediction of the dipole model is 
scaled with an additional free parameter Bl 

Fl{x, Q") = Ff^'\x, Q") (1 + Bt) . (39) 

With Bl as a formal free parameter the GBW fit returns Bt deviating from by more than 3 
standard deviations, Bi = 0.54 + 0.15. The fit for the EM model does not yield a significant 
change for the F^ prediction: 5^ = 0.15 + 0.14. 
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To summarise, a steeper rise to smaller x of the structure function F2, together with a larger 
R value as obtained in the fractal model fit with constant i?, gives the best description of the HI 
data. However, a sufficiently softer rise of Fj together with a smaller Fl, as predicted by the 
IIM model, also describes the data well. For the GBW model, the rise of F2 is rather steep such 
that the fit to the data prefers a larger Fi, which is inconsistent with the prediction of the model. 



12 Summary 

A new measurement is performed of the inclusive double differential cross section for neutral 
current deep inelastic positron-proton scattering, e'^p e^X, in the region of small Bjorken x 
and low absolute momentum transfers squared, Q^. The data were obtained with the HI detector 
at the ep collider HERA in two dedicated periods of data taking at beam energies = 27.5 GeV 
and Bp = 920 GeV. In the year 1999, events were collected with a dedicated trigger on low 
DIS events at the nominal interaction vertex position, corresponding to an integrated luminosity 
of 2.1 pb ^ In the year 2000, the interaction vertex was shifted forward by +70 cm in proton 
beam direction to access even smaller values of Q^, and data with an integrated luminosity of 
0.505 pb"' was taken. 

The measurement is performed in a wide range of inelasticity y, from 0.0015 to 0.8, and 
of Bjorken X, from 5 ■ 10^^ to 0.02. The data cover a range from 0.2 to 12GeV^, with an 
overlap region of the nominal and the shifted vertex data of 0.5 < 2^ < 3.5 GeV^, in which both 
measurements agree. At low the data analysed here comprise the full statistics collected with 
the HI experiment at 920 GeV. 

The measurement obtained with the 1999 and the 2000 data is combined with data col- 
lected in the years 1995 and 1997, which were taken at 820 GeV proton beam energy in similar 
experimental conditions and published previously. This combination takes the correlation of 
systematic uncertainties into account and provides a new, single data set from the HI exper- 
iment, which supersedes all HI data previously released in that kinematic region. The total 
uncertainty of the final reduced cross section measurement is about 2% for a large part of the 
phase space. 

The neutral current ep cross section at low is governed by two independent proton struc- 
ture functions, F2 and F^. For y < 0.6, the influence of the longitudinal structure function Fi 
is small, and the data in this range are also presented as a measurement of the proton structure 
function Fjix, Q^). For y = 0.735, using a method based on the derivative of the cross section 
with respect to In y, the structure function Fl(x, Q^) is extracted with minimum assumptions on 
the behaviour of Fj. 

In each bin a simple parameterisation of the reduced cross section in terms of a power 
law of Fjix, Q^) oc x~'^ and R = Fl/{F2 - Fl) describes the data well. The power A increases 
approximately logarithmically with at > 2GeV^. The parameterisation is consistent 
with a constant value of R{x, Q^) - 0.5, which implies that Fi{x, Q^) - F2{x, Q^)/3 under the 
assumption of a power law rise of F2 towards low x. 

The transition region of DIS to photoproduction, - I GeV^, cannot be analysed within 
perturbative QCD. The data therefore are studied here within phenomenological models. The 



43 



structure function F2{x, Q^) is analysed using a self similarity based ansatz within a fractal 
model. The fractal parameterisation, combined with a constant R, provides a good descrip- 
tion of the measured cross section in the full range of phase space covered by the data. 

The Colour Dipole Model predicts both structure functions F2 and F^ using a single charac- 
teristic dipole scattering cross section. Two versions of the CDM, the GBW model and the IIM 
model, are used in this analysis and are found to generally describe the cross section data well. 
The description of the data in the GBW model is observed to improve when the contribution 
of Fi within this model is formally allowed to be enhanced. The IIM model prediction on Fi 
is similar to the GBW model. Owing to a softer rise of F2 towards small x, the IIM ansatz 
yet is able to describe the cross section data better and no modification on the predicted Fl is 
suggested by the data. 

For the region - I GeV^, in which the transition from photoproduction to DIS takes 
place, the data as presented in this paper are the most precise result of the HI Collaboration. 
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A Averaging Procedure 



The function of equation 28 is to be minimised with respect to the sets m' and bj. This 
determines the averaged measurements and uncertainties, pi''^"^. A, ave, ff/ave and the matrix A^, 
used in equation 29. 

The minimum;^^jj^ in equation 28 is found by solving a system of linear equations obtained 
by requiring dx^/dm' = and dx^/dbj = which can be written in matrix form 



Am Asm \ / ] ^ Cm 
(AsmV As '[ 5-^= I 1 C 



s 



(40) 



Here the vector M^™ corresponds to all measurements and the vector 5^™ corresponds to all 
systematic error sources. The matrix has a diagonal structure with Nm diagonal elements 

e i,e 

The other matrices have the following elements 



"■SM ' 






i,e 

k,e 


C' = 


- y Jh- • 

Zje *2 "^',6' 


d -- 





(42) 



Here 5,;, is the standard Kronecker symbol. Note that the matrix As m has the dimension Nm x Ns 
while the matrix As is quadratic with Ns x Ns elements. 

Using the method of the Schur complement, the solution is found as 

A'^ = As - (AsmVA-^Asm 

5-^= = iA'^r'(Cs-(AsMyA-JCM) (43) 
M-^ = AjjiCM-AsMB^n. 

Given the components of the vector B'^^^, fi^ave = Q'/ave/Ao,^, the solution for /j.'^'"^ can be written 
in explicit form 

fi'^^^^ = ^. (44) 



The uncorrelated uncertainty squared is determined by the inverse of the elements of the diag- 
onal matrix Am 

A^ = — - — (45) 

i.ave 7(1. „ V^-'J 



1 

A 



e a2 



i,e 
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Similarly, the contributions from statistical and systematical uncertainties can be calculated 



a2 _ a 4 ^^\2 a2 _ a4 !^\2 (Af^\ 

!, ave, Stat (,ave / ; ^4 (,e,stat ' r,ave,unc (,ave / ; ^4 i,e,unc" V^*-*/ 



i,e e i.e 



Equation 44 and equation 45 reproduce the standard formula for a statistically weighted average 
of several uncorrelated measurements when all shifts of the systematic error sources are set to 
zero. 

The non-diagonal nature of the matrix expresses the fact that the original sources of the 
systematic uncertainties are correlated with each other after averaging. The matrix can be 
decomposed to re-express equation 27 in terms of diagonalised systematic error sources 

DD = UA's Ya,e = AsmA-JD-' . (47) 

Here U is an orthogonal matrix composed of the eigenvectors of A^ , D is a diagonal matrix 
with corresponding square roots of eigenvalues as diagonal elements and Ya„e represents the 
sensitivity of the average result to these new sources. Its elements are the F:^^*^ in equation 30. 
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0.55 


-1.91 


1.43 


-1.05 


-0.73 


-0.04 


1.20 


1.757 X 


10" 


-5 


6.726 X 


10" 


-1 


0.652 


5.82 


2.66 


2.51 


-1.08 


-0.35 


1.33 


-0.26 


2.16 


-3.57 


1.20 


2.000 X 


10" 


-5 


5.907 X 


10" 


-1 


0.686 


3.95 


2.59 


2.51 


-0.73 


-0.46 


0.40 


-0.04 


0.91 


-0.93 


1.20 


3.200 X 


10" 


-5 


3.692 X 


10" 


-I 


0.697 


3.78 


1.66 


2.73 


-0.81 


-1.73 


0.22 


-0.10 


0.09 


-0.59 


1.20 


6.310 X 


10" 


-5 


1.872 X 


10" 


-1 


0.653 


4.12 


1.37 


2.71 


-1.17 


-2.50 


0.07 


-0.22 


0.22 


-0.02 


1.20 


1.580 X 


10" 


-4 


7.478 X 


10" 


-2 


0.498 


4.40 


2.06 


2.78 
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Table 10: Reduced cross section cr^, as measured with the SVX data sample for 0.2 < < 
1.2 GeV^. The uncertainties are quoted in % relative to cr^.. ^tot is the total uncertainty de- 
termined as the quadratic sum of systematic and statistical uncertainties, ^stat is the statistical 
uncertainty, ^uncor represents the uncorrelated systematic uncertainty. Je'^, jg^, y^^^^, Tnoise. 7e'' 

SpaCal 

and are the bin-to-bin correlated systematic uncertainties in the cross section measurement 
due to uncertainties in the SpaCal electromagnetic energy scale, electron scattering angle, LAr 
calorimeter hadronic energy scale, LAr calorimeter noise, SpaCal hadronic energy scale and the 
photoproduction background, respectively. The luminosity uncertainty of 3% for the SVX data 
is not included in ^tot- 
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Table 11: Reduced cross section cr,., as measured with the SVX data sample for 1.5 < 2^ < 
3.5 GeV^. The uncertainties are quoted in % relative to cr,.. 5tot is the total uncertainty de- 
termined as the quadratic sum of systematic and statistical uncertainties. 5stat is the statistical 
uncertainty, ^uncor represents the uncorrelated systematic uncertainty, y^^, jg^, y£,,^,j, Tnoise, 7e'' 
and jyp are the bin-to-bin correlated systematic uncertainties in the cross section measurement 
due to uncertainties in the SpaCal electromagnetic energy scale, electron scattering angle, LAr 
calorimeter hadronic energy scale, LAr calorimeter noise, SpaCal hadronic energy scale and the 
photoproduction background, respectively. The luminosity uncertainty of 3% for the SVX data 
is not included in ^tot- 
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Table 12: Reduced cross section cr,., as measured with the NVX-BST data sample for 0.5 < 
< 3.5 GeV^. The uncertainties are quoted in % relative to cr,.. ^tot is the total uncertainty 
determined as the quadratic sum of systematic and statistical uncertainties, ^stat is the statistical 
uncertainty, ^uncor represents the uncorrected systematic uncertainty, y^^, jq^, ys,,^^, ynoise> Je'' 
and jyp are the bin-to-bin correlated systematic uncertainties in the cross section measurement 
due to uncertainties in the SpaCal electromagnetic energy scale, electron scattering angle, LAr 
calorimeter hadronic energy scale, LAr calorimeter noise, SpaCal hadronic energy scale and the 
photoproduction background, respectively. The luminosity uncertainty of 1.1% for the NVX 
data is not included in ^w- 
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Table 13: Reduced cross section cr,-, as measured with the NVX-BST data sample for 5 < 
< 12 GeV^. The uncertainties are quoted in % relative to cr^. 5tot is the total uncertainty 
determined as the quadratic sum of systematic and statistical uncertainties, ^stat is the statistical 
uncertainty, ^uncor represents the uncorrelated systematic uncertainty, je'^, je,, y^had' Tnoise, 7e'' 

SpaCal 

and jyp are the bin-to-bin correlated systematic uncertainties in the cross section measurement 
due to uncertainties in the SpaCal electromagnetic energy scale, electron scattering angle, LAr 
calorimeter hadronic energy scale, LAr calorimeter noise, SpaCal hadronic energy scale and the 
photoproduction background, respectively. The luminosity uncertainty of 1.1% for the NVX 
data is not included in Otot- 



54 





X 




\l 


(Ty 






*^uncor 


ye. 


ya. 




Tnoise 


SpiiCal 
% 


yyp 


GeV^ 










% 


% 


% 


% 


% 


% 


% 


% 


1.50 


1.853 X 10" 


-5 


8.000 X 10"' 


0.605 


12.0 


3.18 


3.48 


0.49 


0.81 


-0.02 


0.20 


1.19 


-10.97 


2.00 


2.470 X 10" 


-5 


8.000 X 10"' 


0.756 


9.23 


2.35 


2.70 


-1.47 


1.36 


-0.03 


0.12 


2.15 


-7.98 


2.50 


3.088 X 10" 


-5 


8.000 X 10"' 


0.837 


7.11 


2.46 


2.67 


-1.17 


-0.44 


-0.05 


0.21 


0.12 


-5.98 


3.50 


4.323 X 10" 


-5 


8.000 X 10"' 


0.871 


7.99 


3.10 


2.83 


-0.86 


-0.56 


0.62 


-0.04 


2.32 


-6.28 


5.00 


6.176 X 10" 


-5 


8.000 X 10"' 


0.993 


7.70 


3.12 


2.78 


-1.70 


-0.72 


-0.72 


0.14 


0.50 


-6.13 


6.50 


8.029 X 10" 


-5 


8.000 X 10"' 


1.080 


6.42 


3.11 


2.64 


-0.55 


0.91 


0.00 


-0.46 


-1.62 


-4.54 


8.50 


1.050 X 10" 


-4 


8.000 X 10"' 


1.174 


6.22 


3.73 


2.80 


-0.28 


1.06 


-0.89 


0.00 


-1.12 


-3.71 



Table 14: Reduced cross section a,-, as measured with the NVX-S9 data sample. The uncer- 
tainties are quoted in % relative to cr,.. S^^i is the total uncertainty determined as the quadratic 
sum of systematic and statistical uncertainties, ^stat is the statistical uncertainty, ^uncor represents 
the uncorrelated systematic uncertainty, y^;, je^, Jeu-^^, Tnoise, 7e'' ^^'^ Yyp bin-to-bin 

SpaCiil 

correlated systematic uncertainties in the cross section measurement due to uncertainties in 
the SpaCal electromagnetic energy scale, electron scattering angle, LAr calorimeter hadronic 
energy scale, LAr calorimeter noise, SpaCal hadronic energy scale and the photoproduction 
background, respectively. The luminosity uncertainty of 1.1% for the NVX data is not included 
in Stat. 
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0.100 X 


10 




0.084 


36 


0.85 


0.251 X 


10 




0.033 


37 


0.85 


0.800 X 


10" 




0.010 


38 


0.85 


0.320 X 


10" 




0.003 


39 


1.2 


0.176 X 


10" 




0.675 


40 


1.2 


0.200 X 


10" 


* 


0.593 


41 


1.2 


0.200 X 


10 




0.663 


42 


1.2 


0.320 X 


10" 




0.371 


43 


1.2 


0.320 X 


10" 




0.414 


44 


1.2 


0.631 X 


10" 




0.188 


45 


1.2 


0.800 X 


10" 




0.148 


46 


1.2 


0.130 X 


10" 




0.091 


47 


1.2 


0.158 X 


10 




0.075 


48 


1.2 


0.398 X 


10 




0.030 


49 


1.2 


0.130 X 


10" 




0.009 


50 


1.2 


0.500 X 


10" 




0.002 


51 


1.5 


0.185 X 


10" 




0.800 


52 


1.5 


0.220 X 


10" 




0.675 


53 


1.5 


0.320 X 


10" 




0.463 


54 


1.5 


0.320 X 


10" 




0.518 


55 


1.5 


0.500 X 


10 




0.296 


56 


1.5 


0.800 X 


10 




0.185 


57 


1.5 


0.130 X 


10 




0.114 


58 


1.5 


0.200 X 


10 




0.074 


59 


1.5 


0.320 X 


10 




0.046 


60 


1.5 


0.500 X 


10 




0.030 


61 


1.5 


0.800 X 


10 




0.019 


62 


1.5 


0.100 X 


10 


2 


0.015 


63 


1.5 


0.320 X 


10 




0.005 


64 


1.5 


0.130 X 


10 


I 


0.001 



0.08 
0.06 
0.09 
0.07 
0.06 



0.12 
0.10 
0.08 
0.07 



0.16 
0.15 
0.13 
0.11 
0.10 
0.08 



0.20 
0.19 
0.17 
0.15 
0.13 
0.11 
0.09 



0.22 
0.22 
0.20 
0.20 
0.18 
0.15 
0.13 
0.11 



0.26 
0.26 
0.23 
0.22 
0.21 
0.20 
0.17 
0.14 
0.12 



0.29 
0.29 
0.27 
0.25 
0.23 
0.22 
0.20 
0.19 
0.17 
0.17 
0.14 
0.11 



0.230 
0.190 
0.300 
0.191 
0.203 
0.450 
0.357 
0410 
0.264 
0.237 
0.204 
0.449 
0.442 
0.461 
0.478 
0411 
0.296 
0.280 
0.183 
0.479 
0.504 
0.466 
0.510 
0.678 
0.500 
0.376 
0.308 
0.225 
0.565 
0.614 
0.612 
0.596 
0.567 
0.546 
0.499 
0414 
0.350 
0.307 
0.608 
0.671 
0.741 
0.689 
0.705 
0.647 
0.594 
0.543 
0.503 
0.502 
0.374 
0.298 
0.610 
0.702 
0.756 
0.801 
0.759 
0.699 
0.643 
0.615 
0.584 
0.548 
0.495 
0.463 
0.409 
0.327 



0.230 
0.190 
0.300 
0.191 
0.203 



0.411 
0.264 
0.237 
0.204 



0.472 
0.480 
0.411 
0.296 
0.280 
0.183 



0.490 
0.538 
0.681 
0.500 
0.376 
0.308 
0.225 



0.641 
0.634 
0.573 
0.549 
0.500 
0.414 
0.350 
0.307 



0.714 
0.738 
0.652 
0.597 
0.544 
0.504 
0.502 
0.374 
0.298 



0.804 
0.864 
0.775 
0.705 
0.644 
0.616 
0.584 
0.548 
0.495 
0.463 
0.409 
0.327 



14.3 
13.1 
9.84 
10.00 
10.8 
21.7 
5.74 
9.12 
9.62 
8.81 
9.93 
5.42 
3.75 
19.0 
10.1 
8.85 
8.37 
5.92 
114 
3.96 
3.74 
3.09 
3.02 
17.5 
5.14 
6.79 
4.94 
5.81 
2.54 
5.20 
1.96 
2.65 
1.65 
2.92 
2.78 
5.88 
4.61 
4.56 
2.54 
2.62 
3.60 
1.67 
2.68 
1.18 
2.18 
2.43 
1.67 
2.88 
3.58 
4.51 
3.17 
1.94 
1.77 
1.20 
1.06 
1.26 
1.49 
2.40 
1.60 
2.51 
2.35 
5.22 
2.32 
3.99 



ave.uncor 




CME 


% 


% 


GeV 


12.0 


19.98 


319 


6.18 


15.03 


319 


11.3 


16.82 


319 


4.70 


12.05 


319 


5.29 


12.37 


301 


12.8 


25.34 


319 


11.0 


13.50 


301 


11.1 


20.36 


319 


4.38 


10.99 


319 


4.19 


10.08 


319 


4.55 


11.08 


319 


5.74 


9.44 


319 


9.17 


10.69 


301 


9.84 


21.61 


319 


6.07 


16.25 


319 


4.87 


10.57 


319 


4.20 


9.74 


319 


3.44 


7.07 


319 


6.39 


13.12 


301 


2.90 


5.85 


319 


8.21 


9.89 


301 


5.44 


6.51 


319 


7.28 


8.33 


301 


11.2 


21.16 


319 


5.84 


10.70 


319 


3.46 


7.98 


319 


3.02 


6.17 


319 


3.15 


6.76 


319 


2.52 


4.50 


319 


9.45 


12.17 


301 


5.36 


5.99 


319 


4.98 


6.27 


301 


3.39 


4.13 


319 


4.52 


5.97 


319 


3.59 


5.98 


319 


2.98 


7.31 


319 


2.66 


5.60 


319 


2.81 


5.49 


301 


2.14 


4.65 


319 


2.51 


3.94 


319 


8.36 


9.98 


301 


2.73 


3.55 


319 


4.55 


5.83 


301 


2.25 


3.09 


319 


4.02 


5.24 


319 


4.97 


5.78 


319 


2.30 


3.24 


319 


2.67 


4.26 


319 


2.62 


4.74 


319 


2.60 


5.47 


319 


3.48 


7.93 


319 


1.78 


3.31 


319 


2.12 


3.08 


319 


3.20 


4.47 


301 


1.97 


2.62 


319 


2.15 


2.95 


319 


2.42 


3.32 


319 


2.59 


3.97 


319 


2.18 


3.30 


319 


7.05 


7.74 


319 


2.47 


3.80 


319 


3.74 


6.61 


319 


2.03 


3.51 


301 


2.49 


7.00 


319 



Table 15: Combined HI reduced cross section crf^'^ for 0.2 < 2^ < 1.5 GeV^. The uncertainties 
are quoted in % relative to crf^^. FJ^ represents the structure function Fl used for the CME 
correction (equation 32) and to calculate the structure function F2. 5ave,stat ((5ave,uncor) represents 
the statistical (uncorrelated systematic) uncertainty. 5ave,tot is the total uncertainty calculated as a 
sum of uncorrelated uncertainty and all correlated sources in quadrature. A global normalisation 
uncertainty of 0.5% is not included in ^ave.tot- CME stands for the centre-of-mass energy of the 
measurement. 
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# 

T 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 



yi yi yi y^ y^ y^ y'l n yh ym y\\ y\i y\i yu yis ytt y\i y\% y^ yio ''21 yn y^ y^ yis 

%c %,. %,. %,. 'U %r "/a %r %r %o 

-847 4.05 13.89 -47.47 28.59 -29.31 841 9.60 11.37 0.70 2.45 10.13 4.13 4.53 -7.00 16.08 5.55 -2.12 5.17 -4.79 -1.39 0.57 3.25 -0.76 -1.67 

-8.44 0.03 2.05 -16.70 4.06 3.02 8.85 4.87 -8.92 4.13 7.94 1.45 17.70 1.42 -13.72 -13.88 -0.90 1.39 -7.65 2.47 -9.60 -0.56 1.05 -1.06 -1.12 

-8.46 5.18 15.49 -60.99 27.90 -21.88 11.74 5.01 3.26 2.87 4.43 8.68 9.12 2.88 -10.26 4.44 1.89 -3.94 2.73 -2.79 -3.82 0.68 2.74 -0.82 -1.69 

-8.45 -1.19 1.31 -12.83 11.91 -7.49 20.51 8.35 -7.96 -2.77 22.68 -11.42 11.97 -7.76 -6.84 -1241 0.96 2.99 -8.42 3.25 -10.50 -1.33 2.29 -1.41 -1.41 

-8.46 -1.54 -0.42 -4.36 -2.37 8.37 2.08 9.28 -6.33 2.73 6.68 -1.39 10.58 -1.36 -7.56 -11.91 0.62 1.23 -5.17 1.79 -4.26 -0.14 0.76 -0.69 -0.87 

-8.45 1.47 4.81 -19.08 3.08 1.98 -1.34 3.89 0.63 4.12 -1.29 2.26 -4.14 -1.06 -1.34 -0.65 1.81 -2.91 2.88 -1.95 1.69 0.89 1.76 -0.02 -1.03 

-8.66 19.34 9.32 4.74 8.68 741 -10.14 -4.02 -0.07 0.70 26.17 10.09 -11.30 25.55 5.77 -5.90 1.03 -1.40 -1.20 1.14 0.98 6.38 7.17 -15.10 2.23 

-8.42 7.01 28.05 -123.68 41.80 -15.78 17.04 18.13 -7.89 -1.98 22.54 -2.54 22.13 -12.02 -7.29 -21.45 -2.31 -11.93 4.21 0.70 4.74 2.76 0.80 -0.93 -1.77 

-8.46 0.28 4.59 -17.69 9.30 -6.18 5.69 10.17 0.15 1.78 6.25 0.08 -4.02 -1.68 -5.16 7.06 1.47 0.40 0.05 -2.06 -7.68 -1.15 3.32 -0.83 -1.23 

-8.46 0.26 0.04 -5.57 2.86 0.54 11.44 -1.48 -6.32 3.04 8.08 -3.79 4.03 -2.11 -641 -8.73 0.29 1.19 -5.28 1.42 -8.00 -0.65 2.12 -0.71 -1.22 

-8.47 -0.34 -0.56 -0.84 0.85 1.46 7.59 1.31 -3.73 2.29 6.72 -4.19 -0.09 -2.91 -3.31 -5.60 141 1.01 -3.33 0.59 -5.04 -0.41 2.10 -0.52 -1.12 

-8.44 0.29 9.73 -35.80 6.26 3.70 -10.82 23.30 2.89 2.91 -0.79 4.32 -13.10 -3.38 -0.91 13.08 0.79 -4.09 7.84 -5.04 0.10 -0.19 2.90 -0.32 -0.84 

-8.66 21.09 3.43 2.40 0.05 1.81 -4.02 -3.78 0.51 0.33 20.31 9.38 -7.66 15.14 2.38 -7.80 4.12 -1.96 -0.15 -6.68 1.36 -7.56 1.17 -8.39 1.49 

-8.48 3.46 0.47 -0.95 2.62 -4.76 11.15 -19.63 1.55 3.96 -0.89 -4.72 -11.40 -3.11 5.69 -5.85 3.61 -2.83 1.84 -0.76 4.21 1.70 2.20 0.56 -1.39 

-8.43 4.57 21.52 -96.02 29.98 -7.24 11.43 18.94 -7.85 -1.25 19.35 -4.28 13.28 -11.98 -4.93 -16.82 -2.15 -941 3.50 0.34 2.78 1.81 1.20 -0.81 -1.49 

-8.47 -0.86 3.21 -14.11 7.55 -4.49 3.90 13.42 1.00 -1.32 9.85 -1.87 13.02 -2.64 -3.47 -8.49 441 0.24 -1.94 0.79 2.78 1.05 0.52 -0.60 -1.05 

-8.47 0.90 1.81 -10.87 6.19 -4.16 10.70 -1.34 -3.16 2.24 4.90 -0.34 7.40 1.23 -4.87 -5.96 -0.04 7.70 -5.29 -1.76 -8.17 -0.62 5.40 0.89 0.46 

-8.46 -0.12 -0.41 -1.73 0.12 2.30 6.67 0.72 -6.06 6.10 441 -0.39 -2.26 -0.93 -6.12 -0.98 -0.77 4.02 -2.33 -0.09 -3.58 -0.06 0.94 0.06 0.31 

-8.48 -0.13 0.14 1.36 1.68 -3.36 1.87 2.48 2.66 2.81 0.02 0.58 -0.46 0.58 -1.66 -441 1.11 0.75 2.51 -0.03 -0.55 -0.28 5.49 3.61 4.55 

-8.45 0.84 5.86 -22.11 3.12 3.62 -4.22 10.18 1.17 3.16 -0.05 0.34 -11.86 -4.00 1.27 4.83 1.37 -3.47 5.19 -3.16 1.36 0.38 2.47 -0.04 -0.96 

-8.65 20.32 3.93 2.78 2.24 3.05 -5.19 -3.74 0.38 0.55 20.73 9.06 -8.13 16.77 2.97 -7.18 3.44 -1.80 -0.34 -4.92 1.25 -4.46 2.44 -9.73 1.63 

-8.47 2.20 1.96 -7.04 1.11 -0.54 2.14 -6.90 1.55 5.53 -4.53 3.52 -2.39 2.13 -1.85 -2.62 2.79 -1.88 1.19 -1.35 1.67 1.22 1.57 0.16 -1.13 

-8.60 17.09 -8.57 -1.73 -3.98 -1.64 0.24 -1.41 1.38 2.84 5.70 2.20 -3.48 4.08 -0.37 -3.64 3.89 -0.68 0.22 -2.30 0.72 0.31 4.07 -11.34 1.74 

-8.47 -1.45 -0.74 1.62 0.81 2.74 8.57 5.51 -1.63 -4.27 16.47 -20.04 -17.65 -17.26 13.57 -7.04 348 -1.72 2.01 0.29 5.31 0.52 2.48 0.08 -1.08 

-8.47 2.45 6.38 -31.57 23.73 -22.45 25.69 -2.50 -1.25 -4.12 20.69 -6.10 33.65 -3.25 -6.44 -26.34 5.80 -0.20 -7.40 4.79 4.15 2.58 -0.51 -0.84 -1.72 

-8.47 0.69 1.95 -10.69 5.12 -2.52 8.31 1.20 -2.17 2.61 5.27 -1.29 6.70 0.25 -4.87 -7.55 3.59 4.91 -3.17 -0.93 3.90 1.35 3.44 1.90 1.61 

-8.47 -0.53 -0.73 0.23 1.70 0.53 9.59 2.56 -4.82 2.94 6.83 -4.26 -7.21 -2.28 -7.14 -0.15 -1.88 3.65 -1.02 0.16 -6.45 -0.93 2.24 0.34 0.86 

-8.47 -0.35 -0.56 0.36 -0.63 2.74 4.23 1.95 -1.82 0.93 3.68 -2.31 -2.55 -1.35 -4.59 -2.15 0.31 1.73 1.17 -0.82 5.09 1.38 0.84 1.12 1.28 

-8.45 1.53 5.16 -18.70 3.31 1.29 -2.22 4.47 1.99 3.96 -1.96 1.47 -11.17 -2.20 0.92 4.49 1.87 -3.25 4.75 -3.07 1.44 0.59 2.45 0.05 -1.04 

-8.70 23.84 10.36 5.53 6.16 5.84 -9.16 -5.70 0.01 -1.55 30.75 15.38 -12.02 23.47 4.49 -9.88 3.79 -2.20 -0.78 -6.07 1.57 -6.03 2.46 -10.56 1.81 

-8.47 2.26 1.04 -3.27 0.34 -0.80 3.24 -9.51 1.66 548 -4.58 2.61 -3.13 1.92 -0.93 -3.44 3.08 -1.71 0.93 -1.14 2.10 1.32 1.56 0.24 -1.15 

-8.61 18.05 -7.99 -1.01 -1.26 -0.16 -0.97 -1.35 0.94 2.93 7.86 1.83 -3.62 7.39 0.73 -378 3.22 -1.06 0.16 -2.86 0.79 -1.16 3.13 -9.74 1.53 

-8.48 4.42 -1.99 -0.46 -0.97 -0.59 2.40 -5.61 1.56 3.93 -1.05 1.99 -2.33 0.01 -1.63 -3.62 4.04 -0.48 0.25 -0.76 1.14 2.56 3.85 -7.64 0.64 

-8.47 -0.30 2.10 0.44 -3.54 -1.65 3.26 -4.98 241 1.35 -0.36 6.63 -3.37 -9.21 -5.03 -4.90 7.60 0.55 0.67 -3.16 0.67 0.65 5.53 -15.11 2.29 

-8.47 0.31 4.11 -5.19 -6.12 1.27 1.30 5.71 -0.09 12.30 -5.40 -16.02 -0.58 17.29 6.25 5.63 -3.93 3.08 -2.00 9.65 -0.45 16.99 7.93 -6.56 1.04 

-8.47 0.70 1.64 -10.06 6.23 -4.71 9.66 0.29 -1.57 1.96 5.60 0.03 16.26 2.56 -7.16 -13.29 4.58 5.94 -6.66 1.56 2.20 1.15 1.81 1.22 1.08 

-8.47 -0.07 0.12 -2.35 1.87 -0.19 7.08 1.88 -3.63 4.38 4.05 -0.74 -3.24 0.11 -6.71 0.99 0.65 4.88 -3.20 043 -4.44 -0.60 1.12 0.15 0.54 

-8.47 0.36 0.10 1.05 1.75 -3.55 2.87 -0.56 2.01 2.88 -0.51 2.68 0.72 2.00 -3.79 -0.34 1.29 0.86 0.55 0.04 -0.93 0.21 1.36 0.75 0.95 

-8.45 2.26 6.77 -24.90 4.39 1.50 -2.50 3.92 241 4.89 -3.89 1.71 -12.77 -0.87 0.71 4.36 4.00 -0.73 1.74 -1.45 1.19 -0.02 1.97 0.46 -0.33 

-8.46 0.65 1.86 -6.42 -0.75 3.25 -1.89 2.34 1.30 4.00 -2.04 1.28 -5.72 -0.73 -0.22 -0.07 2.54 -1.65 2.11 -1.76 1.69 0.75 1.74 0.05 -0.96 

-8.65 19.02 3.66 3.59 5.94 4.61 -6.00 -4.30 0.16 0.33 20.90 9.48 -8.11 16.48 2.88 -7.48 3.63 -1.84 -0.29 -5.55 1.28 -5.67 1.86 -8.97 1.54 

-8.47 1.68 1.25 -4.31 0.72 -0.48 2.54 -5.52 1.46 4.70 -2.90 2.12 -2.03 1.26 -1.32 -3.13 3.06 -1.39 0.71 -1.06 1.71 1.13 1.53 0.11 -1.11 

-8.58 13.51 -6.83 -1.18 -1.40 -0.20 -0.45 -1.46 1.28 3.32 4.69 1.21 -3.48 4.63 -0.05 -2.66 3.21 -0.40 0.09 -0.24 0.55 3.91 5.50 -12.79 1.88 

-8.47 2.22 0.22 0.06 -1.15 -1.11 4.10 -9.48 1.56 4.51 -3.05 2.71 -3.06 -0.46 -1.70 -4.29 4.16 -0.91 041 -1.56 1.40 0.86 2.51 -3.84 -0.21 

-8.47 -1.72 1.85 0.72 1.37 1.19 0.80 -3.89 1.92 2.81 -0.58 3.26 -3.55 -3.91 -2.92 -2.39 5.00 0.80 0.22 1.74 0.32 8.73 8.37 -17.40 2.47 

-8.47 -0.17 1.27 -0.73 -1.94 -0.14 4.34 0.00 0.50 6.35 -2.71 -8.02 2.23 3.73 048 -0.47 1.76 -1.37 0.86 -3.27 0.54 -3.55 0.52 -4.62 0.76 

-8.47 0.91 0.45 -0.64 -0.55 0.22 2.09 -3.37 0.86 6.94 -5.11 2.27 -0.47 4.82 -3.90 -0.38 2.49 1.82 -0.81 -1.70 -1.85 0.01 2.70 -1.26 0.15 

-8.47 -0.39 -0.73 -1.57 3.31 1.37 5.23 1.37 -1.41 1.38 4.91 -6.20 4.31 -0.08 -0.36 -6.73 1.58 1.83 -1.94 0.69 1.77 1.75 2.49 -1.38 0.83 

-8.47 -0.04 -0.37 0.82 2.94 -3.09 8.60 0.29 -1.85 2.71 5.70 -2.84 -4.57 -1.95 -3.82 0.16 0.37 3.44 -0.22 -0.06 -0.56 0.34 0.56 0.06 0.21 

-8.48 0.32 0.61 2.45 3.93 -7.98 2.06 1.24 8.00 -1.04 -0.84 0.21 -0.57 0.75 1.25 -3.12 0.95 2.57 441 -0.82 2.32 0.86 1.35 0.67 0.48 

-8.45 3.77 9.23 -29.52 3.30 2.81 -5.20 2.53 11.83 5.11 -16.00 -3.97 -18.62 7.34 10.88 -6.95 16.43 35.69 -19.63 5.70 0.80 -3.19 -0.64 1.91 1.75 

-8.46 1.44 3.59 -12.55 1.39 1.58 -1.29 1.02 2.21 4.21 -3.19 0.65 -8.34 -0.10 1.71 1.08 4.70 3.06 -1.91 0.49 1.05 -0.20 0.58 0.42 0.05 

-8.46 0.80 1.49 -4.84 -0.54 2.39 -0.68 0.56 2.13 2.79 -2.49 0.96 -4.78 0.54 -1.40 -0.62 3.37 -0.89 0.46 -0.71 1.37 0.50 1.37 0.17 -0.67 

-8.65 20.32 12.14 4.13 1.88 6.94 0.93 -2.78 -0.14 3.49 3.77 -1.86 2.24 -1.59 -2.25 0.39 6.11 1.73 4.71 -0.14 -1.24 0.27 1.25 -0.11 -2.26 

-8.47 4.02 0.50 -0.36 -0.62 0.76 2.62 -5.38 -0.07 3.58 -0.77 0.76 -0.68 0.72 -1.28 -2.75 4.33 -0.12 1.68 -0.61 0.73 0.90 1.79 -1.70 -1.00 

-8.47 0.86 0.50 -0.07 -1.58 0.22 3.02 -7.42 2.14 4.55 -2.01 2.48 -1.88 -0.90 -2.34 -3.74 4.70 -0.50 0.81 -1.75 0.88 0.83 2.14 -3.35 -043 

-8.47 1.76 -0.40 048 1.50 0.05 4.06 -9.57 0.86 5.55 -3.82 1.09 -1.97 2.02 -148 -3.47 3.18 -1.28 0.34 -2.79 0.70 -2.24 0.56 -0.14 -0.77 

-8.47 0.33 0.67 0.60 -1.50 -1.74 6.77 -2.79 0.96 2.27 1.91 -5.61 -4.44 -5.07 2.10 -3.96 4.39 -1.39 1.25 -4.56 2.18 -5.80 -0.34 -0.16 -0.53 

-8.47 1.15 0.52 -2.31 0.17 -1.24 6.11 -6.12 -0.41 4.87 -0.93 0.56 6.18 1.02 -3.14 -9.45 3.90 0.73 -1.76 -1.99 2.02 -2.19 0.74 0.40 0.18 

-8.47 -0.02 -0.71 0.26 1.13 0.47 4.47 -1.74 0.49 4.78 -2.20 -4.58 2.27 -0.83 -1.30 -2.88 4.07 -1.58 1.25 -7.91 0.86 -11.42 -2.45 -1.80 0.50 

-8.47 0.52 -0.26 0.19 2.39 -2.63 8.11 -3.04 -0.48 2.09 4.58 -3.30 -0.35 -1.97 -1.10 -5.22 3.09 0.36 -1.91 0.19 -1.19 0.42 1.85 -0.24 -1.19 

-8.47 -0.04 -0.31 0.03 0.90 -0.21 5.50 0.66 -3.43 5.31 3.31 0.36 -3.80 0.06 -7.25 1.77 1.88 0.65 -0.95 1.72 2.21 0.26 0.44 1.87 3.16 

-8.47 -0.03 -0.22 2.39 3.48 -5.32 6.97 1.21 0.52 3.37 4.01 -1.22 -5.44 -1.03 -3.25 2.64 1.71 1.73 0.51 -0.30 -2.30 -0.11 1.42 0.01 -0.02 

-8.50 0.05 2.81 4.38 3.81 -14.06 -15.68 9.14 25.44 -5.63 -8.31 3.97 9.46 0.23 18.29 -8.77 11.73 -4.45 10.74 -0.87 22.17 3.93 -0.06 1.60 0.11 



e 16: Correlated systematic uncertainties given in per mill relative to the cross section measurement reported in table 15. The coefficients 
72™ represent diagonalised correlated systematic uncertainties (see section 9.1). 
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Table 17: Combined HI reduced cross section crj'^'^ for 2 < 2^ < 5 GeV^. The uncertainties 
are quoted in % relative to cr^^'^ . F^l^ represents the structure function Fl used for the CME 
correction (equation 32) and to calculate the structure function F2. 5ave,stat ((^ave.uncor) represents 
the statistical (uncorrelated systematic) uncertainty. 5ave,tot is the total uncertainty calculated as a 
sum of uncorrelated uncertainty and all correlated sources in quadrature. A global normalisation 
uncertainty of 0.5% is not included in 5ave,tot- CME stands for the centre-of-mass energy of the 
measurement. 
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nm 
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-8.47 
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nil 
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-8.48 
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0.68 


7 QO 


..46 


7 13 

-7.43 


0.68 


^ 1 o 

5.19 


4.35 


3.65 


n 07 


7 


3 nft 
-3.06 


1 3n 
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3 77 
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0.14 
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-1.61 


0.05 
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0.10 


n 3 
-0. . 


n 73 
0.7. 
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0.1! 


0.5; 


n 1/1 
0.44 


1.53 


105 


-8.46 


1 .80 


4.92 


-14.97 


2.59 


-1.94 


-2.92 


2.30 


1.72 


!4.68 


-5.70 


1.97 


-9.50 


0.82 


9.43 


-2.10 


1 ! .20 


17.94 


-11.17 


5.22 


0.25 


-1.82 


-1.08 


1.70 


2.67 


-0.99 


106 


-8.46 


0.50 


1.24 


-4.03 


0.10 


0.05 


-0.67 


-0.67 


0.33 


3.67 


-0.4! 


1.63 


-2.4! 


-0.96 


2.76 


1.49 


2.55 


3.52 


-3.39 


2.73 


0.06 


-0.40 


-2.26 


0.00 


0.83 


-0.17 


107 


-8.66 


19.06 


14.15 


5.00 


3.27 


3.75 


-3.04 


-2.46 


0.80 


3.95 


1.40 


-5.04 


1.88 


-5.00 


-3.91 


-0.38 


7.23 


0.85 


2.97 


-2.86 


-2.95 


4.10 


-0.70 


1.67 


-2.10 


-1.62 


108 


-8.46 


0.04 


0.14 


-0.41 


-0.13 


-0.54 


0.28 


-0.74 


-0.67 


2.96 


1.41 


2.03 


-0.0! 


-1.63 


1.27 


1.54 


0.48 


-0.21 


-0.63 


1.62 


-0.07 


-0.01 


-1.03 


0.39 


!.29 


-0.09 


109 


-8.51 


2.68 


2.06 


0.99 


0.37 


-0.45 


-0.48 


-0.85 


0.05 


2.70 


1.99 


-0.99 


-0.54 


-1.07 


-1.83 


-1.01 


3.54 


0.06 


1.14 


-1.47 


-1.32 


1.98 


0.36 


1.11 


-2.17 


-1.38 


110 


-8.46 


-0.04 


0.00 


0.10 


-0.23 


-0.48 


0.02 


-0.94 


1.92 


3.66 


2.58 


0.32 


-0.60 


-1.15 


-1.10 


-0.31 


1.76 


-0.26 


0.36 


-0.29 


-0.74 


1.08 


0.20 


0.90 


-1.19 


-1.02 


Hi 


-8.46 


-0.05 


-0.0! 


0.6! 


0.35 


-1.54 


0.83 


0.12 


-1.36 


2.69 


0.93 


1.92 


-0.43 


-1.04 


-0.39 


0.65 


-0.41 


-0.39 


-0.19 


2.04 


0.67 


-1.07 


1.36 


0.55 


-1.14 


-0.75 


112 


-8.47 


-0.17 


-0.13 


0.42 


0.10 


-1.25 


0.56 


-0.32 


-0.92 


3.65 


1.67 


1.70 


-0.43 


-1.17 


-0.75 


0.83 


0.54 


-0.74 


-0.12 


2.12 


0.51 


-0.91 


1.3! 


0.70 


-1.10 


-0.84 


113 


-8.47 


-0.20 


-0.40 


0.66 


-0.05 


-0.45 


1.66 


-0.58 


-0.76 


2.03 


1.70 


1.40 


0.58 


-0.79 


-1.28 


0.20 


-1.27 


0.84 


0.77 


0.67 


0.92 


-0.84 


3.05 


1.21 


-0.81 


-0.89 


114 


-8.47 


-O.iO 


-0.57 


0.25 


-0.22 


1.30 


3.40 


-0.59 


0.98 


3.05 


2.16 


2.04 


0.21 


-0.91 


-1.49 


0.95 


-1.13 


0.84 


0.89 


2.09 


l.!6 


-1.42 


2.54 


0.88 


-1.46 


-0.95 


115 


-8.47 


-0.30 


-0.47 


-0.23 


-0.72 


0.43 


1.45 


-1.12 


-2.73 


0.58 


1.64 


0.76 


0.18 


0.15 


-0.83 


0.34 


2.86 


0.57 


2.39 


1.24 


-0.30 


-0.08 


1.12 


0.90 


-1.02 


-0.89 


116 


-8.47 


-0.07 


-0.26 


1 .33 


1.09 


-1.28 


3.77 


1.55 


0.41 


4.00 


1.27 


2.10 


-1.76 


-0.70 


-2.08 


-0.30 


-1.38 


0.41 


0.59 


1.64 


0.85 


-0.61 


1.06 


0.65 


-1.22 


-0.90 


117 


-8.48 


-0.12 


0.86 


2.97 


3.00 


-7.66 


-1.74 


4.35 


5.83 


4.88 


- 1 .40 


1.72 


-3.17 


-0.97 


0.24 


-4.68 


-1.40 


-1.78 


1.56 


2.23 


-1.66 


-0.25 


-0.32 


0.67 


0.69 


-0.50 



Table 18: Correlated systematic uncertainties given in per mill relative to the cross section measurement reported in table 17. The coefficients 
yT^' 72™ represent diagonalised correlated systematic uncertainties (see section 9.1). 
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iSavc.slal 


iSavc.iiiiior 


(Savc.tol 


CME 




GeV^ 












% 


% 


% 


GeV 


118 


6.5 


0.803 X 10 


0.800 




1.083 




3.11 


2.64 


4.85 


319 


119 


6.5 


0.951 X lO""* 


0.675 




1.053 




2.95 


2.31 


4.16 


319 


120 


6.5 


0.130 X 10"^ 


0.494 


0.45 


1.123 


1.211 


1.67 


2.22 


2.97 


319 


121 


6.5 


0.130 X 10"-* 


0.552 


0.45 


1.124 


1.239 


1.53 


1.73 


3.03 


301 


122 


6.5 


0.200 X 10"-* 


0.321 


0.41 


1.123 


1.152 


1.25 


2.20 


2.72 


319 


123 


6.5 


0.200 X 10"-^ 


0.359 


0.41 


1.117 


1.155 


1.09 


1.62 


2.24 


301 


124 


6.5 


0.320 X 10"^ 


0.201 


0.37 


1.006 


1.015 


0.84 


1.16 


1.75 


319 


125 


6.5 


0.500 X 10"^ 


0.128 


0.34 


0.936 


0.939 


0.86 


1.26 


1.86 


319 


126 


6.5 


0.800 X 10"-* 


0.080 


0.31 


0.854 


0.855 


0.87 


1.26 


1.91 


319 


127 


6.5 


0.130 X 10"- 


0.049 


0.28 


0.758 


0.758 


0.90 


1.28 


1.87 


319 


128 


6.5 


0.200 X 10-- 


0.032 


0.26 


0.694 


0.694 


0.92 


1.29 


1.89 


319 


129 


6.5 


0.398 X 10-- 


0.016 


0.22 


0.616 


0.617 


0.69 


1.19 


1.76 


319 


130 


6.5 


0.130 X 10-' 


0.005 


0.18 


0.482 


0.482 


0.73 


1.80 


2.44 


319 


131 


8.5 


0.105 X 10-3 


0.800 




1.178 




3.72 


2.80 


5.20 


319 


132 


8.5 


0.124 X 10-' 


0.675 




1.211 




2.26 


2.28 


3.44 


319 


133 


8.5 


0.139 X 10-' 


0.675 




1.136 




2.07 


1.82 


4.54 


301 


134 


8.5 


0.200 X 10-' 


0.420 


0.46 


1.178 


1.239 


1.52 


2.22 


2.88 


319 


135 


8.5 


0.200 X 10-' 


0.469 


0.46 


1.182 


1.261 


1.38 


1.64 


2.59 


301 


136 


8.5 


0.320 X 10-' 


0.262 


0.41 


1.112 


1.131 


0.91 


1.25 


1.86 


319 


137 


8.5 


0.500 X 10-' 


0.168 


0.37 


1.033 


1.039 


0.95 


1.18 


1.81 


319 


138 


8.5 


0.800 X 10-' 


0.105 


0.34 


0.950 


0.953 


0.95 


1.28 


1.90 


319 


139 


8.5 


0.130 X 10-- 


0.065 


0.30 


0.842 


0.842 


0.99 


1.30 


1.94 


319 


140 


8.5 


0.200 X 10-- 


0.042 


0.28 


0.773 


0.773 


1.00 


1.30 


1.93 


319 


141 


8.5 


0.320 X 10-- 


0.026 


0.25 


0.663 


0.663 


1.04 


1.32 


1.99 


319 


142 


8.5 


0.631 X 10-- 


0.013 


0.22 


0.604 


0.604 


0.79 


1.24 


1.83 


319 


143 


8.5 


0.200 X 10-' 


0.004 


0.17 


0.456 


0.456 


0.88 


1.82 


2.67 


319 


144 


12.0 


0.800 X 10-' 


0.148 


0.38 


1.053 


1.058 


1.07 


1.30 


1.99 


319 


145 


12.0 


0.130 X 10-- 


0.091 


0.34 


0.923 


0.924 


1.10 


1.30 


1.97 


319 


146 


12.0 


0.200 X 10-- 


0.059 


0.31 


0.861 


0.861 


1.11 


1.33 


2.00 


319 


147 


12.0 


0.320 X 10-- 


0.037 


0.28 


0.757 


0.757 


1.14 


1.34 


2.02 


319 


148 


12.0 


0.631 X 10-- 


0.019 


0.24 


0.646 


0.646 


0.88 


1.24 


1.86 


319 


149 


12.0 


0.200 X 10-' 


0.006 


0.19 


0.490 


0.490 


0.93 


1.83 


2.51 


319 



Table 19: Combined HI reduced cross section cr^^'^ for 6.5 < <\2 GeV^. The uncertainties 
are quoted in % relative to af"^. represents the structure function Fi used for the CME 
correction (equation 32) and to calculate the structure function F2. ^ave.stat ((5ave,uncor) represents 
the statistical (uncorrelated systematic) uncertainty. 5ave,tot is the total uncertainty calculated as a 
sum of uncorrelated uncertainty and all correlated sources in quadrature. A global normalisation 
uncertainty of 0.5% is not included in 5ave,tot- CME stands for the centre-of-mass energy of the 
measurement. 
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'I 
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'3 
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%c 
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'6 

%c 


V 


'8 
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%c 


ni 

%c 


9 
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%c 


'15 

%p 


'16 
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%c 
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%o 
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%o 


yVf 

'23 

%o 
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'24 


'25 

%o 
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118 


-8.46 


1 37 


3 17 


-1 1 27 


1 39 


1 70 


1 22 


86 


2 94 


3 14 


-4 40 


-1 73 


-7 26 


3 19 


65 


-2 69 


1 29 


14 18 


-7 90 


1 22 


1 47 


-1.01 


3.45 


3 17 


3.62 


-1.11 


119 


-8 46 


93 


2 26 


-7 40 


54 


40 


-1 09 


-0 24 


1 49 


4 29 


-2 28 


96 


-4 56 


08 


3 74 


48 


4 45 


7 88 


-5 60 


3 07 


26 


-0.74 


-1.96 


33 


1 05 


-0.37 


120 


-8 46 


1 1 


30 


-1 03 


-0 17 


-0 41 


00 


-0 89 


-0 86 


3 43 


1 27 


2 20 


-0 25 


-1 76 


1 91 


1 85 


70 


43 


-1 18 


1 85 


-0 04 


-0 04 


-1 59 


15 


1 05 


-0 06 


121 


-8 59 


10 47 


7 68 


2 70 


1 52 


2 48 


-1 22 


-1 95 


1 85 


4 01 


2 36 


-2 94 


75 


-3 10 


-2 97 


-0 64 


5 55 


49 


2 18 


-2 22 


-2 23 


3 14 


-0 14 


1 46 


-2 26 


-1 59 


122 


-8 46 


01 


08 


-0 12 


-0 06 


-0 69 


35 


-0 60 


-0 53 


2 91 


1 46 


1 98 


-0 03 


-1 59 


96 


1 40 


56 


-0 69 


-0 43 


1 7 1 


-0 1 1 


-0 03 


-0 89 


47 


1 39 


-Oil 


123 


-8 50 


1 34 


1 04 


58 


03 


-0 48 


-0 36 


-1 00 


99 


3 18 


2 42 


-0 64 


-0 63 


-0 98 


-1 87 


-0 91 


2 88 


06 


93 


-1 14 


-1 05 


1 60 


65 


1 07 


-2 30 


- 1 39 


124 


-8 46 


-0 45 


-0 72 


18 


-0 13 


-1 37 


-0 10 


-0 54 


-0 01 


2 96 


2 05 


49 


-0 69 


-0 83 


-0 57 


-0 43 


2 06 


-0 36 


38 


-0 49 


-0 84 


1 24 


-0 05 


83 


-0 97 


-0 93 


125 


-8 47 


42 


34 


90 


62 


-1 60 


89 


28 


-0 78 


3 33 


85 


2 24 


-0 29 


-1 48 


-0 21 


78 


-1 42 


-0 30 


-0 55 


2 22 


93 


-1 36 


1 57 


48 


- 1 24 


-0 72 


126 


-8 47 


-0 16 


-0 12 


55 


17 


-1 64 


40 


-0 31 


-7 32 


4 68 


1 84 


2 28 


-0 07 


-1 54 


-0 07 


1 30 


68 


-0 91 


-0 21 


2 51 


64 


-1 14 


1 36 


68 


-0 96 


-0 79 


127 


-8.47 


-0.15 


-0.39 


0.67 


0.16 


-0.34 


2.57 


-0.18 


-1.07 


1.93 


1.56 


1.46 


-0.25 


-0.56 


-1.76 


0.28 


-1.50 


1.45 


1.07 


0.88 


0.99 


-0.63 


1.90 


0.81 


-1.09 


-0.82 


128 


-8.47 


-0.20 


-0.55 


-0.10 


-0.69 


0.95 


2.17 


-1.19 


-2.02 


0.46 


1.55 


1.40 


0.82 


-0.21 


-0.21 


0.48 


0.35 


1.92 


2.57 


1.03 


0.19 


-0.46 


143 


0.64 


-1.41 


-0.79 


129 


-8.47 


-0.03 


-0.42 


0.92 


0.64 


0.12 


4.65 


1.14 


-0.15 


2.55 


1.11 


2.21 


-1.06 


-0.39 


-1.29 


-0.07 


-0.97 


0.94 


1.13 


1.29 


0.67 


-0.46 


1.06 


0.59 


-1.42 


-0.89 


130 


-8.48 


-0.19 


0.52 


2.59 


2.43 


-6.37 


-0.78 


3.29 


4.07 


4.47 


-0.50 


1.58 


-2.72 


-0.82 


-0.58 


-3.79 


-2.07 


0.11 


2.60 


1.96 


-1.09 


-0.17 


-0.54 


0.50 


0.44 


-0.44 


131 


-8.47 


1.14 


2.93 


-8.35 


2.27 


-0.86 


0.79 


2.55 


5.45 


2.64 


-4.75 


-1.59 


-7.97 


3.23 


-0.57 


-4.39 


5.86 


11.44 


-4.64 


2.62 


0.47 


-1.14 


1.96 


2.49 


3.18 


-0.93 


132 


-8.46 


0.46 


1.11 


-3.62 


0.02 


0.09 


-0.64 


-0.79 


0.37 


3.16 


-0.27 


1.63 


-2.11 


-1.03 


2.63 


1.55 


2.05 


3.14 


-2.99 


2.60 


0.02 


-0.33 


-2.38 


-0.12 


0.69 


-0.12 


133 


-8.70 


24.85 


18.62 


7.45 


5.55 


2.79 


-3.21 


-0.87 


-2.55 


2.64 


-1.56 


-4.85 


2.24 


-6.11 


-3.80 


0.10 


5.30 


0.87 


2.28 


-1.70 


-2.06 


2.95 


-0.09 


141 


-2.01 


-1.37 


134 


-8.46 


0.06 


0.15 


-0.56 


-0.20 


-0.36 


0.25 


-0.85 


-0.81 


3.03 


1.48 


2.02 


-0.12 


-1.67 


1.20 


1.85 


0.85 


-0.55 


-0.99 


1.84 


-0.02 


-0.05 


-1.19 


0.34 


1.25 


-0.09 


135 


-8.54 


6.02 


4.54 


1.91 


1.12 


0.06 


-0.36 


-0.51 


-2.38 


1.40 


1.01 


-1.83 


-0.24 


-1.17 


-1.65 


-1.24 


5.60 


0.13 


1.95 


-2.34 


-2.10 


3.03 


-0.41 


1.22 


-1.85 


-1.37 


136 


-8.46 


-0.08 


0.01 


0.06 


-0.27 


-0.79 


-0.14 


-0.85 


0.14 


2.98 


2.32 


0.20 


-0.56 


-0.89 


-0.58 


-0.36 


2.91 


-0.49 


0.52 


-0.78 


-1.14 


1.65 


-0.36 


0.92 


-0.76 


-0.94 


137 


-8.47 


-0.47 


-0.26 


0.05 


-0.28 


-0.98 


-0.11 


-0.77 


0.43 


2.90 


2.22 


0.42 


-0.56 


-0.87 


-0.56 


-0.38 


2.06 


-0.29 


0.42 


-0.59 


-0.89 


1.31 


-0.08 


0.84 


-0.97 


-0.93 


138 


-8.46 


-0.09 


-0.04 


0.61 


0.28 


-1.64 


0.72 


0.02 


-2.19 


3.17 


1.14 


2.06 


-0.30 


-1.12 


-0.16 


0.86 
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Table 20: Correlated systematic uncertainties given in per mill relative to the cross section measurement reported in table 19. The coefficients 
yT^' represent diagonalised correlated systematic uncertainties (see section 9.1). 
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Table 21: Orthogonal transition matrix Uj^ from the original (b^', etc) to the diagonalised (Z?j™,Z72™) systematic sources for the averaged HI 
data. The matrix ellements are given in %. 




Figure 1: Illustration of the bins used for the cross section measurement in the NVX (left) and 
SVX (right) analyses. Dashed lines of constant 9g indicate approximate angular acceptance for 
both measurements. The dark (light) shaded area corresponds to the bins where the electron (S) 
method is used for the measurement of the cross section. The measurement in the bins outside 
the angular acceptance range employ the 2 method for ISR events. 





Figure 2: A low event as reconstructed in the HI detector. The electron is scattered into 
the backward region. The electron trajectory is reconstructed in the Backward Silicon Tracker 
(BST) and in the Backward Drift Chamber (BDC). The electron energy is determined using 
the SpaCal calorimeter. The hadronic final state is detected in the central and forward tracking 
detectors, and in the LAr calorimeter. 
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Figure 3: Schematic layout of the HI Backward Silicon Tracker (BST). The active area is 
composed of eight wheels subdivided into two modules, BSTl and BST2, of four wheels each. 
One wheel is made of 16 r sensors and one u sensor (mounted on the back side, not shown 
here). Eight consecutive sensors in 1/16 of azimuth build a BST sector. In z the module BSTl 
extends from -73.2 to -95.7 cm, BST2 from -35.9 to -46.9 cm. Readout boards are placed in 
the rear section. Also indicated are the electric shielding and the water cooling pipes. 




Figure 4: The two types of silicon sensors used in the BST: a) r sensor, b) u sensor, each with 
640 readout strips. The r sensor has a double metal structure for the readout lines to reach the 
top (outer radius) part where the five amplifiers are mounted on the hybrid, as sketched. 
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Figure 5: Three-dimensional view of a section of the BDC illustrating the layer structure and 
the drift cell geometry. The chamber has a radial coverage from 6 cm to 71 cm. At a radius of 
about 22 cm the segmentation is changed and a transition drift cell is introduced. 
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Figure 6: Efficiency of the SpaCal electron triggers, S9 a), SO b), and S3 c), used in this analysis, 



as a function of E'^. 
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Figure 7: Central Tracker vertex reconstruction efficiency as a function of y-^. The dashed line 
corresponds to the applied selection criterion, > 0.03. In the analysis the cross section at 
high y is measured with y^ instead of y^- 
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Figure 8: a) BST track segment finding efficiency as a function of the radial position of the 
electron candidate in the SpaCal, for the NVX data sample, b) ratio of data to MC efficiencies. 
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Figure 9: a) BST track segment finding efficiency as a function of tlie radiai position of tlie 
eiectron candidate in tlie SpaCal, for tlie SVX sample, b) ratio of data to MC efficiencies. 
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Figure 10: Relative alignment of the Central Tracker (CT) and BDC. Left: ^ct - ^bdc versus 
^CT after alignment. Right: ^ct - 6'bdc after alignment. 
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Figure 11: a) Scatter plot of the number of BST hits linked to the electron candidate BST 
track as a function of the azimuthal angle 0e determined by the associated SpaCal cluster. At 
least three linked hits are required to define a track. A number of linked hits exceeding eight 
corresponds to a track passing the azimuthal BST wafer overlap region; b) Contours of equal 
density for the distribution of A9 = 9\ - 62, where 0i 2 are the polar angles measured in the two 
overlapping BST sectors, as a function of 0e- The horizontal dotted lines indicate +0.2 mrad as 
is used for the systematic uncertainty of the alignment. 
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Figure 12: a) Distribution of the scattered electron energy E'^ for the NVX data sample; 
b) Kolmogorov-Smirnov test probability distribution as a function of the relative shift in the 
measured and simulated energy distributions. 




Figure 13: a) Energy distribution for n candidates based on the NVX sample as triggered by the 
low energy trigger, S9; b) di-photon invariant mass distribution for tjP candidates. The double 
angle calibration constants are applied to the data and MC simulation. 
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Figure 14: Invariant mass distribution of the two electron candidate tracks for a special J/ifr 
event selection. The line indicates a fit to the data. Mfn and o-fn correspond to the Gaussian 
mean and width of the peak. 
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Figure 15: Summary of SpaCal energy scale determination. The band indicates the uncertainty 
due to the scale diff'erence between the data and the simulation. 
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Figure 16: Distributions of the scattered electron energy E[, for the data and the MC simulation 
in the NVX (left) and the SVX (right) analyses. The MC bands include the statistical uncertainty 
and the effect of a +0.2% electromagnetic energy scale variation. 
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Figure 17: Transverse momentum balance Pj/Pj distribution for the data and the MC simu- 
lation in the NVX and the SVX analyses. The bands include the statistical uncertainty of the 
simulation and the effect of the LAr hadronic scale uncertainty, see description in the text. The 
vertical line indicates the analysis requirement P^/P^ > 0.3. 
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Figure 18: E-P- distribution for the data and the MC simulation in the NVX and the SVX 
analyses. The bands include the statistical uncertainty of the simulation and the effect of a 
+0.5 GeV variation of the SpaCal hadronic final state contribution. 
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Figure 19: Relative contributions to the measured from the LAr (closed circles), tracks (trian- 
gles) and SpaCal (open circles) together with the subtracted LAr noise fractions (squares) in the 
NVX (left) and SVX (right) analyses. The distributions of simulated events are shown as curves. 
The shaded areas correspond to a 10% systematic uncertainty on the LAr noise description. 
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Figure 20: Distribution of E'^ a) and 6e b) for photoproduction events detected in the electron 
tagger. The plots are based on the NVX-S9 analysis. The dashed line in a) corresponds to the 
minimum E'^ permitted by the analysis cut i/e < 0.85. 
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Figure 21: Distribution of E'^ for background events, estimated using wrong charge BST tracks 
(equation 25) for data and the PHOJET simulation. The simulated sample is normalised using 
photoproduction events with the scattered electron detected in the electron tagger. 
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Figure 22: Distribution of events for the NVX-BST analysis: the energy a) and the polar angle 
b) of the scattered positron; E-P^ c) and the z vertex position d). Control distributions for 
the NVX-S9 analysis: energy e) and polar angle f) of the scattered positron candidates. The 
histograms represent the simulation of DIS and the photoproduction background (shaded). 
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Figure 23: Distribution of events for the SVX-BST a-c,e,f) and SVX-BDC d) analyses: the 
energy a) and the polar angle b) of the scattered positron; E-P- c) and the z vertex position d); 
the energy e) and E-P^ f) for the ISR bins. The histograms represent the simulation of DIS and 
the photoproduction background (shaded). 
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Figure 24: Distribution of Bjorken-x and using the electron (top) and sigma (bottom) recon- 
struction methods for the NVX data. The histograms represent the simulation of DIS and the 
photoproduction background (shaded). 
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Figure 25: Distribution of Bjorken-x and using the electron (top) and sigma (bottom) recon- 
struction methods for the SVX data. The histograms represent the simulation of DIS and the 
photoproduction background (shaded). 
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Figure 26: Comparison of reduced cross sections as obtained with the electron (closed circles) 
and S (open circles) reconstruction methods, for the NVX data sample. The errors represent 
statistical uncertainties only. 
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Figure 27: Reduced inclusive e^p scattering cross section as measured in the NVX-BST (open 
circles), NVX-S9 (triangles) and SVX (closed circles) analyses of the 920 GeV data. The errors 
represent the statistical and systematic uncertainties added in quadrature. 
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Figure 28: Reduced cross section cTr- Closed circles: combined 1999-2000 data taken at Ep = 
920 GeV; Triangles: SVX data taken in 1995 [35]; Open circles: NVX data taken in 1997 [37]. 
The normalisation of the 1997 data has changed by +3.4%, see section 8.4. The 1995 and 1997 
data were taken at Ep = 820 GeV but are corrected here for comparison to 920 GeV. The errors 
represent the statistical and systematic uncertainties added in quadrature. 
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Figure 29: Reduced cross section cr^, from the combined low HI data, as a function of x 
compared to the X fit result (solid line) and to a /I parameterisation with the same values of 
c{Q^) and A{Q^) but R = (dashed line). The errors represent the statistical and systematic 
uncertainties added in quadrature. 
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Figure 30: Reduced cross section cr^-, from the combined low HI data, as a function of 
X compared to the GBW and IIM models. The errors represent the statistical and systematic 
uncertainties added in quadrature. 
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Figure 31: Structure function F2, from the combined low HI data for y < 0.6, as a function 
of X compared to the fractal, the dipole GBW and the dipole IIM fit results. The errors represent 
the statistical and systematic uncertainties added in quadrature. 
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Figure 32: Measurement of the virtual photon-proton cross section cr^f^ as a function of 
at various values of W. The cross sections for different W values are multiplied with the fac- 
tors indicated in the figure. The errors represent the statistical and systematic errors added in 
quadrature. The averaged HI results are compared to data obtained by the ZEUS experiment 
and to the fractal and dipole model fit results. 
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Figure 33: Coefficients c and A, as defined in equation 34, determined from a fit to the HI 
data as a function of Q^. The inner error bars represent uncorrelated systematic uncertainties. 
The outer error bars represent total uncertainties. The line in b) shows a straight line fit for 
Q^>2 GeV2. 
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Figure 34: Coefficient i? as a function of from a simple parameterisation of the reduced cross 
section as defined inequation 34. The dashed line is drawn at = 0.5. The errors represent the 
total uncertainties. 
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Figure 35: Derivative i/avAcr^/Ai/ for the combined 1999-2000 HI data compared to the predic- 
tions of the dipole models and the fractal model for F2 with an assumption R = 0.5 to describe 
Fl, labeled R = 0.5. The lines increasing as a function of Iny correspond to Ft = for these 
models. The lines turning over at high y correspond to the cross section predictions. The in- 
ner error bars represent statistical and uncorrelated uncertainties added in quadrature, the outer 
error bars represent the total uncertainties. 
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Figure 36: Structure function Fi extracted using the derivative method. The solid line is drawn 
for R = 0.5 assuming the fractal parameterisation for Fj. The dashed (dotted) line corresponds 
to the dipole GBW (IIM) model. The inner error bars represent statistical and uncorrected 
uncertainties added in quadrature, the outer error bars represent the total uncertainties. The 
solid (yellow) band indicates the model uncertainty, see text. 
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Figure 37: Comparison of the structure functions F2 (left) and Fi (right) for = 1.2 GeV^ as 
a function of Bjorken jc, for the fractal fit with R = 0.5 (solid line), and the predictions of the 
dipole models, GBW (dashed line) and IIM (dotted line), resulting from the fits to the HI cross 
section data. The vertical line indicates the value of x = Xg for which the GBW dipole model 
saturation radius is reached. 
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